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Executive  Summary 

This  report  reviews  current  approaches  and  provides  adjusted  methodologies  for  predicting 
the  probability  of  blast-related  injuries  following  nuclear  detonations  in  both  urban  and 
open  field  environments.  The  purpose  of  this  research  is  to  integrate  the  optimal  models 
into  our  Health  Effects  of  Nuclear  and  Radiological  Environments  platform  which  is  a  suite 
of  physiological  models  that  predicts  health  effects  in  nuclear  and  radiological  environments. 
The  blast  injury  criteria  models  predict  the  likelihood  of  tertiary  and  secondary  injuries 
based  on  blast  wave  and  object  parameters.  Separate  models  are  used  for  different  injury 
levels,  including  moderate,  serious,  and  fatal  categories. 

Current  approaches  for  predicting  the  likelihood  of  blast-related  injury  use  a  two-step 
process.  First,  the  maximum  velocity  obtained  by  the  human  body  or  secondary  missile 
is  determined  through  the  interpolation  of  values  in  look-up  tables.  These  tables  provide 
velocity  as  a  function  of  peak  overpressure  and  yield.  Next,  probit  models  are  used  to 
predict  the  probability  of  secondary  or  tertiary  injury  as  a  function  of  the  object’s  velocity. 
We  sought  to  review  both  these  steps  and,  if  needed,  adjust  the  methodologies  to  predict 
injury  in  urban  environments.  The  traditional  method  to  determine  an  object’s  velocity 
used  by  the  nuclear  effects  communities  assumes  a  specific  blast  wave  shape.  However, 
in  urban  environments,  complex  overpressure  and  dynamic  pressure  waveforms  are  likely. 
Thus,  rather  than  using  the  peak  overpressure  and  yield  to  calculate  velocity,  we  propose 
the  use  of  a  translation  model  that  uses  physical  principles  to  relate  an  object’s  acceleration 
to  the  dynamic  pressure  and  wind  velocity  as  a  function  of  time.  For  the  next  step  of  the 
process  (i.e.,  relating  velocity  to  injury),  we  review  existing  probit  models  used  by  nuclear 
effects  communities.  This  review  includes  a  detailed  analysis  of  source  data,  assumptions, 
and  ambiguities  in  each  probit  model  generation.  In  several  cases  where  the  models  were 
based  on  limited  data,  included  unrealistic  assumptions,  and/or  large  ambiguities  existed, 
we  developed  new  criteria  models  using  contemporary  data.  Finally,  based  on  our  review 
we  recommend  models  to  be  included  in  our  Health  Effects  of  Nuclear  and  Radiological 
Environments  platform. 
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1  Introduction 


Applied  Research  Associates,  Inc.  (ARA)  has  been  tasked  by  the  Defense  Threat  Reduction 
Agency  (DTRA)  to  support  their  mission  to  safeguard  against  weapons  of  mass  destruction. 
A  subtask  of  this  project  involves  implementing  models  that  predict  the  probability  of  mor¬ 
tality  and  injury  following  nuclear  detonations  in  urban  environments.  This  report  reviews 
current  approaches  for  estimating  blast-related  casualties  and,  where  appropriate,  recom¬ 
mends  alternate  methodologies.  The  optimal  approaches  will  be  selected  for  integration  into 
our  Health  Effects  of  Nuclear  and  Radiological  Environments  (HENRE)  platform  which  is  a 
suite  of  physiological  models  that  predict  prompt  and  protracted  health  effects  of  nuclear  and 
radiological  environments.  Ultimately,  HENRE  models  will  be  integrated  into  other  tools 
such  as  NucFast  or  DTRA’s  Hazard  Prediction  and  Assessment  Capability  (HPAC)  tool  to 
improve  current  casualty  estimation  capabilities  and  better  account  for  urban  environments, 
combined  injuries,  and  the  time-course  of  effects.  Both  HPAC  and  NucFast  provide  location 
and  time-specific  blast  wave  parameters  that  can  serve  as  input  to  the  blast  injury  criteria 
models.  Therefore,  based  on  the  reliability  of  the  underlying  data,  we  identify  the  most  suit¬ 
able  blast  injury  criteria  models  for  implementation  into  HENRE  and  subsequent  integration 
into  HPAC  or  NucFast. 

The  three  categories  of  blast  injury  are  primary,  secondary,  and  tertiary.  Primary  injuries 
are  caused  by  the  direct  effects  of  the  blast  wave  overpressure  on  the  human  body.  Secondary 
injuries  are  caused  by  missiles  that  are  accelerated  by  the  blast  wave.  Tertiary  injuries  are 
caused  by  the  acceleration  of  the  human  body  and  the  ensuing  deceleration.  In  this  work, 
we  focus  on  secondary  and  tertiary  injuries.  Because  of  the  dramatic  effects  experienced 
from  blast,  thermal,  and  radiation  interactions  associated  with  nuclear  detonations,  persons 
experiencing  primary  blast  injuries  will  also  encounter  other  lethal  exposures.  Therefore,  for 
our  purposes,  we  exclude  an  analysis  of  primary  blast  injuries. 

Calculating  the  probability  of  secondary  or  tertiary  injury  due  to  blast  is  a  two-step 
process:  (1)  human  or  missile  velocities  are  determined  based  on  blast  wave  parameters 
and  (2)  probit  models  relate  the  velocity  of  an  object  to  the  probability  of  injury.  The 
Personnel  Risk  and  Casualty  Criteria  (PRCC)  report  and  Effects  Manual-1  Chapter  14 
(EM-1  Ch.  14)  provide  starting  points  for  this  endeavor  (Department  of  the  Army  2013; 
Drake  et  ah  1993;  Reeves  2015).  Both  of  these  documents  use  look-up  tables  developed  in 
the  1970s  to  determine  the  velocity  of  a  human  or  missile  based  on  peak  overpressure  and 
yield  (Fletcher  et  ah  1975).  Next,  probit  models  are  used  to  relate  object  velocity  to  the 
likelihood  of  injury. 

Based  on  known  shortcomings  and  a  review  of  the  literature,  we  identified  two  near-term 
issues  that  could  be  addressed  to  more  accurately  predict  blast-related  injuries  in  urban 
environments: 

1.  Velocity  Calculations:  Current  approaches  use  peak  overpressure  and  yield  to  predict 
injury,  but  this  is  not  an  accurate  method  when  considering  the  complex  waveforms 
that  occur  in  urban  environments.  We  suggest  a  revised  approach  in  which  a  trans¬ 
lation  model  developed  by  Lovelace  Foundation  for  Medical  Education  and  Research 
(LFMER)  is  used  to  predict  the  velocity  of  objects  based  on  the  incident-specific  blast 
wave  parameters  (see  Section  3). 
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2.  Injury  Calculations:  The  existing  probits  relating  velocity  to  injury  are  based  on  lim¬ 
ited  data  (Section  4  provides  a  review  of  these  probit  models  and  the  corresponding 
source  data).  Due  to  the  many  assumptions  involved  in  the  generation  of  the  exist¬ 
ing  probit  models,  where  possible,  we  propose  new  probit  models  that  integrate  more 
contemporary  data. 

The  combined  translation  and  probit  models  provide  improvements  to  probability  of  injury 
predictions  for  complex  blast  wave  scenarios. 
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2  Methods 


To  select  models  for  integration  into  HENRE  2.0,  we  reviewed  current  approaches  for  es¬ 
timating  injury  and  mortality  due  to  nuclear  blast.  In  reviewing  existing  methodologies, 
we  focused  on  approaches  currently  used  by  the  nuclear  effects  communities.  It  is  possi¬ 
ble  that  relevant  models  used  by  other  communities  exist;  however,  these  are  not  reviewed 
here.  Sources  reviewed  include  EM-1  Ch.  14  (Drake  et  al.  1993;  Reeves  2015),  the  PRCC 
2013  report  (Department  of  the  Army  2013),  the  Allied  Medical  Publication  8  (AmedP-8) 
North  Atlantic  Treaty  Organization  (NATO)  planning  guide  (NATO  2009),  and  methodolo¬ 
gies  present  by  Dr.  John  Mercier  at  the  Human  Response  Panel  (HRP)  in  the  early  2000s 
(Mercier  2001).  The  models  presented  by  Dr.  Mercier  are  used  within  the  Probability  of 
Damage  Calculator  (PDCALC),  a  tool  used  by  HPAC  to  estimate  human  injury  in  nuclear 
environments  (Jackson  et  al.  2013).  In  reviewing  existing  models,  our  goal  was  to  select 
the  optimal  models  for  inclusion  in  HENRE  2.0.  Thus,  we  sought  to  fully  understand  the 
derivation  and  accuracy  of  existing  models. 

While  reviewing  the  existing  probit  models  used  to  relate  velocity  to  injury,  several 
issues  arose.  In  some  cases,  the  data  used  to  develop  the  models  were  not  obtainable  or  not 
referenced  and,  thus,  could  not  be  reviewed.  In  other  cases,  when  the  data  was  available, 
the  existing  probit  could  not  be  derived.  Finally,  in  several  cases,  potentially  unrealistic 
assumptions  were  used  to  generate  the  probit  models. 

When  these  issues  arose,  we  sought,  if  feasible,  to  generate  new  probit  models  by  either 
re-deriving  the  probit  using  the  original  data  or  creating  a  new  probit  using  contemporary 
data.  In  some  cases,  we  developed  more  general  probit  equations.  For  example,  multiple 
penetration  probits  exist  for  missiles  of  different  masses.  Instead  of  using  these  probits,  we 
developed  probits  which  include  both  mass  and  velocity  in  the  probit  equation. 

All  computational  procedures  were  performed  in  R  v3.1.2  (R  Core  Team  2014).  Trans¬ 
lation  simulations  were  performed  using  the  general  solver  for  ordinary  different  equations 
from  the  deSolve  library  (Soetaert  et  al.  2010b).  Least-square  regression  was  performed  us¬ 
ing  the  modCost  and  modFit  tools  from  the  Flexible  Modeling  Environment  (FME)  Library 
(Soetaert  et  al.  2010a).  Probit  models  were  fit  to  data  using  generalized  linear  models  (Ven¬ 
ables  et  al.  2002).  R  code  for  running  the  translation  model  and  generating  the  output  data 
and  plots  given  in  Section  3  is  provided  in  Appendix  A. 
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3  Translation  Model 


In  the  1960s,  LFMER  developed  a  translation  model  which  predicts  the  acceleration,  velocity 
and  translation  distance  of  objects  exposed  to  a  blast  wave  based  on  attributes  of  the  object, 
including  mass  and  orientation  (Bowen  et  al.  1961).  An  updated  version  of  this  model 
takes  into  account  deceleration  due  to  ground  friction  (Fletcher  et  al.  1966).  Using  these 
models,  LFMER  performed  a  series  of  simulations  to  develop  look-up  tables  for  different 
initial  orientations  which  relate  the  maximum  velocity  of  a  human  to  blast  wave  parameters 
(Fletcher  et  al.  1975).  These  tables  are  referenced  today  by  both  the  2013  PRCC  report 
and  EM-1  Ch.  14  as  the  source  for  blast  injury  criteria  (Drake  et  al.  1993;  Reeves  2015; 
Department  of  the  Army  2013).  Through  the  interpolation  of  the  values  in  these  tables, 
velocity  is  calculated  as  a  function  of  peak  overpressure  and  yield.  The  approach  used  in  the 
2013  PRCC  report  and  EM-1  Ch.  14  assumes  the  blast  waveforms  of  interest  have  the  same 
shape  and  duration  as  the  ones  originally  used  to  generate  the  look-up  tables. 

In  this  section,  we  discuss  an  alternate  approach  for  calculating  an  object’s  velocity. 
To  provide  the  probability  of  injury  after  exposure  to  a  complex  blast  wave,  rather  than 
using  look-up  tables,  we  have  implemented  the  original  translation  model  developed  by 
LFMER  (Bowen  et  al.  1961;  Fletcher  et  al.  1966).  The  maximum  velocity  of  an  object 
is  calculated  using  the  time-dependent  overpressure,  dynamic  pressure,  wind  velocity,  and 
shock  wave  velocity.  This  method  eliminates  assumptions  about  the  shape  and  duration 
of  the  blast  wave  that  are  inherent  when  using  the  look-up  tables.  While  this  approach 
requires  additional  time-dependent  environment  inputs  leading  to  longer  computing  times, 
it  also  improves  accuracy  by  more  precisely  calculating  the  maximum  velocity  for  a  given 
scenario.  This  section  describes  our  implementation  of  the  translation  model,  as  well  as  the 
verification  and  validation  of  the  model. 

In  developing  the  translation  model,  we  assume  the  direction  of  the  dynamic  pressure 
and  wind  velocity  is  constant.  However,  due  to  the  complex  nature  of  urban  environments 
this  may  not  be  the  case.  Later  versions  of  the  model  should  use  vector  notation  to  predict 
how  changes  in  wind  and  dynamic  pressure  direction  would  affect  translation. 
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3.1  Translation  Model  Nomenclature 

In  developing  the  translation  model,  we  use  a  terminology  similar  to  Bowen  et  al.  1961. 
Dimensional  quantities  are  represented  with  a  lowercase  letter  (except  for  temperature).  The 
dimensionless  quantity  is  either  represented  with  the  corresponding  uppercase  character  or, 
if  unavailable,  with  a  tilde.  The  parameters  used  to  make  quantities  dimensionless  do  not 
change  during  the  blast  scenario. 


Time: 

t  =  time  after  blast  wave  arrival 
tp  =  duration  of  positive  overpressure 
t+  =  duration  of  positive  winds 
tp  =  t/t+ 
tu  =  t/t+ 

Velocity: 

Co  =  speed  of  sound  in  undisturbed  air 
u  =  wind  velocity 
x  =  shock  wave  velocity 
v  =  velocity  of  the  moving  object 
U  =  u/cq 
X  =  x/cq 
V  =  v/c0 

Pressure: 

Po  =  ambient  pressure 
p  =  overpressure 
q  =  dyanmic  pressure 
P  =  p/p0 

Q  =  q/Po 


Distance: 

d  =  distance  traveled  by  object 
D  =  d/(cot+) 

Object  Parameters: 

m  =  mass  of  object 

s  =  area  of  object  presented  to  the  wind 
Cd  =  drag  coefficient  of  object 
a  =  sCd/m  (acceleration  coefficient) 
m0  =  average  mass  of  human;  75  kg  (165  lb) 
A  =  ap0t+/c0 
M  =  m/rriQ 

Air  Temperature: 

T  =  temperature 
T0  =  ambient  temperature 
T  =  T/T0 

Air  Density: 

p  =  air  density 
R  =  gas  constant  of  air 
p  =  pRT0/P0 


Note  that  an  s  subscript  refers  to  the  value  of  that  variable  at  the  shock  front.  For  example, 
ps  is  the  overpressure  at  the  shock  front. 
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3.2  Object  Acceleration 

This  section  presents  our  implementation  and  verification  of  the  translation  model  developed 
by  Bowen  et  ah  1961.  This  model  describes  an  object’s  translation  profile  (i.e. ,  acceleration, 
velocity,  and  distance  traveled)  over  time  following  a  blast  wave  exposure.  Appendix  A. 2 
gives  the  R  function  that  was  used  to  run  the  model.  For  small  debris,  this  model  is  sufficient; 
however,  for  larger  objects,  such  as  animals  and  humans,  the  incorporation  of  deceleration 
is  essential  and  will  be  described  in  Section  3.3. 


3.2.1  Model  Implementation 

The  following  equation  gives  the  time-dependent  acceleration  of  a  generic  object  exposed  to 
a  blast  wave  (Bowen  et  al.  1961) 


dv  f  u  —  v  \  2 

-dt=qa{—) 


(1) 


where  v  is  the  object’s  velocity,  q  is  the  dynamic  pressure,  and  u  is  the  wind  velocity.  Note 
that  q  and  u  are  the  values  at  the  location  of  the  object,  and  this  location  changes  as  the 
object  moves  through  space.  The  acceleration  coefficient  a  is  specific  to  the  object  of  interest, 
taking  into  account  attributes  such  as  surface  area  exposed  to  the  blast  wave,  shape,  and 
mass  (see  Section  3.1,  Nomenclature).  This  equation  was  derived  by  equating  the  drag  force 
of  the  blast  wave  winds  with  the  force  needed  to  cause  object  acceleration. 

A  dimensionless  form  of  Equation  1,  given  by  Bowen  et  ah  1961,  can  be  written  as  follows 
with  the  appropriate  substitutions  (see  Section  3.1,  Nomenclature): 


dV 

dtu 


QA 


fU -V 

\U~ 


2 


(2) 


Equation  2  defines  the  acceleration  portion  of  the  translation  model  that  will  be  used  to 
predict  the  maximum  velocity  an  object  obtains. 


3.2.2  Model  Verification 


To  verify  that  the  acceleration  model  (Equation  2)  correctly  matches  the  original  version  of 
the  model  (Bowen  et  ah  1961),  we  used  the  blast  wave  parameters  specified  in  the  report  as 
the  inputs  and  compared  the  resulting  outputs. 

The  time  profile  for  the  dimensionless  dynamic  pressure  at  a  point  in  space  is  defined  as 
(Bowen  et  ah  1961) 

Q  =  Qs(l  ~  tu)(Je~ ^  +  Ke~s*u)  (3) 

where 


Qs 

J 


2.5 Ps2  1  +  2  •  10“8P4 
7  +  PS  1  +  10-®P4 

fl.l86Ps1/3  if  Ps  <  0.6 
1  if  0.6  <PS< 


10  4P71/4 

104+Pf 


if  Ps  >  1.0 


1.0 


(4) 

(5) 
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K  =  1  -  J  (6) 

7  =  1/4  +  3.6  P]/2  (7) 

5  =  7  +  8 Ps1/2  +  2P2  / (240  +  Ps).  (8) 

The  time  profile  for  the  dimensionless  wind  velocity  at  a  point  in  space  is 

U  =  Ua(l-tu)e~^  (9) 

where 

Us  =  Ps/il  +  Pi'2)  (10) 

V  =  Ps1/3  +  0.0032PS3/2.  (11) 


Equations  2-11  define  the  blast  wave  parameters  at  a  constant  point  in  space.  To  run 
model  simulations,  we  are  interested  in  the  blast  wave  parameters  at  the  location  of  the 
object.  To  make  this  transformation  we  change  tu  in  Equations  3  and  9  to 

D 

tu  -  (12) 

where  Xs  represents  the  dimensionless  speed  of  the  shock  wave  propagation  and  D  represents 
the  dimensionless  distance  traveled  by  the  object  at  time  tu.  Equation  12  represents  the 
amount  of  time  since  the  shock  wave  front  was  at  D,  rather  than  the  amount  of  time  since 
the  shock  wave  front  was  at  the  object’s  starting  position.  This  transformation  assumes  that 
any  decay  due  to  radial  expansion  of  the  shock  wave  is  negligible. 

The  dimensionless  velocity  of  the  pressure  propagation  A"  is  related  to  the  dimensionless 
wind  velocity  U  according  to  the  following  Equation1 : 

Xs  =  jj Us  +  (13) 

The  resulting  simulations  agree  with  outputs  tabulated  in  the  report  (Figure  3.1).  For 
all  simulations  shown  in  Figure  3.1,  Ps  =  0.10.  Simulations  were  stopped  when  the  object’s 
velocity  was  equal  to  the  wind  velocity  at  the  location  of  the  object. 

In  a  second  verification  to  dimensional  output  given  by  Bowen  et  al.  1961,  the  parameters 
in  Table  3.1  were  used.  The  relationship  between  £+  and  was  determined  through  the 
digitization  of  Figure  2.2  from  Bowen  et  al.  1961.  Based  on  these  parameters,  a  maximum 
velocity  of  7.01  m/s  (23.0  ft/s)  was  obtained  which  is  very  close  to  the  value  given  in  the 
report  of  7.13  m/s  (23.4  ft/s). 

Based  on  these  verification  results,  we  are  satisfied  that  the  acceleration  model  is  accu¬ 
rately  reproduced  here.  The  slight  discrepancies  between  the  model  and  the  data  that  exist 
are  likely  caused  by  errors  in  digitization  and  different  numerical  procedures. 

1This  equation  was  derived  by  Bowen  et  al.  1961  from  Shapiro  1954  (Pg  1001,  Equation  25.20)  by  setting 

u\  to  zero,  uijc\  to  U,  W/c\  to  X ,  and  k  to  1.4. 


Velocity  (dimensionless) 


Figure  3.1:  Verification  of  the  acceleration  model  implementation.  Data  digitized 
from  Bowen  et  al.  1961,  Table  4.1. 

Table  3.1:  Parameters  used  for  dimensional  acceleration  model  verification. 


Parameter 

Value 

a 

0.0061  m2/kg  (0.03  ft2/lb) 

Ps 

36.5  kPa  (5.3  psi) 

Po 

91.7  kPa  (13.3  psi) 

4 

0.964  s 

ti 

1.216  s* 

Co 

341.4  m/s  (1120  ft/s) 

*Value  was  obtained  through  the  digitiza¬ 
tion  of  Figure  2.2  from  Bowen  et  al.  1961. 
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3.3  Object  Deceleration 

This  section  describes  the  implementation  and  verification  of  the  decelerative  portion  of  the 
translation  model  that  was  developed  by  LFMER  (Fletcher  et  al.  1966).  The  deceleration 
rate  was  quantified  by  LFMER  using  experimental  animal  data  with  initial  speeds  ranging 
from  4.47  to  26.82  m/s.  The  stopping  distance  was  determined  as  a  function  of  the  initial 
velocity  and  the  mass  of  the  animal.  This  empirical  relationship  was  in  turn  used  to  determine 
the  deceleration  rate. 


3.3.1  Model  Implementation 

The  relationship  between  velocity  and  stopping  distance  is 

logi»  (*  © 1 1/i)  =  ~a+i  lo&»  (°  © I/e)  (14) 

where  x  is  the  distance  traveled  before  stopping,  m  is  the  weight  of  the  animal,  mo  is  the 
average  weight  of  a  man  or  75  kg  (165  lbs),  and  v  is  the  initial  velocity.  Using  Equation  14, 
the  instantaneous  deceleration  is  given  as2 


(20) 


where  F  =  4.278  m1"^3-2  and  B  =  0.38308. 

The  deceleration  term  given  by  Equation  20  is  made  dimensionless  using  the  appropriate 
substitutions  (see  Section  3.1,  Nomenclature)  and  is  incorporated  into  the  translation  model 

2  Equation  14  implies  that  the  distance  traveled  over  the  time  interval  [to, tend]  is 

x{tend )  -  a: (t0)  =  10”a  •  (— )  7  7  v{t0)b.  (15) 

We  assume  this  relationship  holds  for  all  t  £  [to, tend]-  That  is  ,  we  assume  the  distance  “remaining”  at  time 
t  (i.e.,  the  distance  traveled  over  the  time  interval  [t,tend])  is 

/  771n  \  1/3 

x{tend)  -  x(t)  =  10-a  •  (  —  )  v(t)b. 

\  nn  / 

Therefore, 

Differentiating  Equation  17  with  respect  to  t,  we  obtain 

and  solve  for  deceleration 

dv  1  (  /  Trio  \  2— k 

dt  b ■ 10_a  \  V  to  /  ) 

Equation  19  provides  the  deceleration  of  an  object  at  instantaneous  velocity  v.  This  result  agrees  with  that 
given  by  Fletcher  et  al.  1966  (a  =  0.8399  log10(mb-1/s)  and  b  =  1.6169). 


(16) 

(17) 

(18) 

(19) 
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given  by  Equation  2.  The  final  equation  for  the  dimensionless  translation  model  is 

dv-cA(u-vy  nj  ( v  y 

3.3.2  Model  Verification 

To  verify  the  final  translation  model  given  by  Equation  21,  we  compare  results  of  our  model 
simulation  with  results  from  a  LFMER  report  (Fletcher  et  al.  1975).  In  this  report  Fletcher 
et  ah  performed  simulations  to  predict  the  translational  profile  of  exposed  personnel.  For 
each  simulation,  the  peak,  positive  duration,  and  impulse  of  both  the  overpressure  and  dy¬ 
namic  pressure  are  given.  Simulations  were  run  for  personnel  in  different  initial  orientations, 
including  prone  and  standing.  Appendix  A. 3  gives  the  R  function  used  to  run  model  simu¬ 
lations  to  perform  the  verifications. 

To  verify  the  model  given  by  Equation  21,  the  dynamic  pressure  wave  and  time-course  of 
the  wind  velocity  must  be  specified.  Assumptions  are  required  to  obtain  these  inputs  because 
the  overpressure  and  dynamic  pressure  waveforms  were  not  made  explicit  in  the  1975  report. 
Thus,  these  waveforms  must  be  estimated  here  using  the  peak,  duration,  and  impulse. 

The  waveforms  for  the  dimensionless  overpressure  and  dynamic  pressure  are  assumed  to 
be  of  the  following  form 


P  =  P8(l  -  tp)e~nt*>  (22) 

Q  =  Qs(l-tu)e~r^.  (23) 


This  is  in  agreement  with  equations  used  in  other  LFMER  reports  to  define  pressure  waves 
(Bowen  et  ah  1968).  Using  Equation  22  and  23,  the  dimensional  overpressure  impulse  (Ip) 
and  the  dynamic  pressure  impulse  (Iq)  for  the  positive  portion  of  the  blast  wave  are 

Ip=—f(e~n  +  n-  1)  (24) 

7lz 

=  +  (25) 

For  a  specific  scenario,  the  value  of  n  and  r  were  determined  by  solving  Equation  24  and  25 
since  the  values  for  ps,  qs,  Ip,  Iq,  tp,  and  are  given. 

To  determine  the  time-course  of  the  wind  velocity,  we  (1)  calculate  the  temperature  at 
the  shock  wave  front  using  Rankine-Hugoniot  relations, 

Ts  =  7  +  Ps  (1  +  PS)  (26) 

s  7  +  6  PSK  J  J 

(2)  calculate  the  temperature  within  the  shock  wave  assuming  adiabatic  conditions  (an  as¬ 
sumption  that  is  stated  in  Fletcher  et  ah  1975), 


T 


Ts 


f  P  +  l\*? 

\Ps  +  l) 


(27) 
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Figure  3.2:  Verification  of  the  final  translation  model.  Model  simulation  (red 
line/circles)  compared  with  digitized  data  (black  line/circles)  obtained  from  Fig¬ 
ure  A-l  of  Fletcher  et  al.  1975.  The  open  circles  connected  by  dashed  lines 
represent  the  same  time  points  post  blast  wave  arrival. 


where  7  is  the  adiabatic  index  for  a  diatomic  gas  or  7/5,  and  (3)  calculate  the  air  density. 

_  P  +  1 
P~  f 

Finally,  the  dimensionless  wind  velocity  can  be  calculated  using  q  =  \pu2  and  the  appropriate 
substitutions  (see  Section  3.1,  Nomenclature): 


U 


2QRTq 

<%P 


(29) 


Equations  22-29  were  used  to  solve  Equation  21  numerically.  Results  comparing  our 
model  output  with  a  simulation  from  Fletcher  et  al.  1975  are  shown  in  Figure  3.2.  This 
simulation  represents  the  10th  run  in  the  Fletcher  1975  report  for  a  prone  personnel  in  a 
random  orientation.  This  verification  represents  a  10  kT  surface  burst  and  a  ground  range 
of  394  m  (1263  ft).  The  parameters  used  for  this  simulation  are  given  in  Table  3.2.  The 
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Table  3.2:  Parameters  used  to  run  verification  shown  in  Figure  3.2. 


Parameter  Value 


Ps 


189  kPa  (27.4  psi) 

0.355  s 

17.1  kPa  s  (2.478  psi  s) 
96.5  kPa  (14  psi) 

0.672  s 

7.05  kPa  s  (1.021  psi  s) 
287.058  J  kg"1  K_1 
298  K 


following  equation  is  used  to  define  a  (as  specified  in  Fletcher  et  al.  1975): 


(30) 


where  d  is  the  distance  the  object  has  traveled  in  meters. 

Figure  3.2  shows  slight  discrepancies  between  our  model  output  and  the  translational 
profile  given  by  Fletcher  et  al.  1975.  Specifically,  the  peak  velocity  predicted  by  our  model 
is  slightly  lower  and  the  resulting  predicted  displacement  is  less.  The  discrepancies  may 
be  caused  by  the  assumed  shape  of  the  dynamic  and  overpressure  waveforms  or  differences 
between  the  numerical  methods.  Given  the  sources  of  error,  we  are  satisfied  that  we  have 
accurately  reproduced  the  translation  model  with  deceleration. 
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3.4  Validation 


We  validated  the  translation  model  by  comparing  two  sets  of  displacement  data  to  model 
predictions.  The  first  set  of  data  relates  the  dynamic  pressure  impulse  Iq  to  the  displacement 
of  jeeps  (Needham  2010),  and  the  second  set  relates  the  displacement  of  pedestrians  who 
have  been  struck  by  automobiles  to  the  speed  of  the  automobile  during  impact  (Otte  2001). 
Appendix  A. 3  gives  the  R  function  used  to  run  model  simulations  to  perform  the  validations. 

3.4.1  Displacement 

The  dynamic  pressure  impulse  is  the  only  blast  wave  parameter  available  for  the  jeep  dis¬ 
placement  data  (Needham  2010).  Therefore,  we  ran  the  model  with  different  blast  wave 
parameters  subject  to  the  constraint  that  the  dynamic  pressure  impulse  had  to  equal  a  spec¬ 
ified  value.  £+  and  were  varied  from  0.1  to  1  s,  r  and  n  were  varied  from  0.1  to  5.  Using 
these  values  and  Equation  4,  22  and  23,  we  determined  Ps  and  Qs.  We  approximated  the 
value  of  a  for  a  jeep  side  on  as  0.005  m2/kg  ( Cd  =  1,  s  =  5  m2,  m  =  1000  kg).  This  process 
was  repeated  for  dynamic  pressure  impulses  ranging  from  0.5  to  100  kPa  s.  Based  on  these 
simulations,  we  determined  the  minimum  and  maximum  displacement  as  a  function  of  the 
dynamic  pressure  impulse. 

Figure  3.3  compares  the  predicted  displacement  ranges  to  the  data  on  jeep  displacement. 
From  this  comparison,  it  is  evident  that  the  model  is  able  to  accurately  predict  jeep  displace- 


Figure  3.3:  Model  prediction  (shaded  region)  of  jeep  displacement  range  as  a 
function  of  dynamic  pressure  impulse.  Data  (red  dots)  from  Needham  2010. 
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ment,  despite  the  many  assumptions  involved  in  the  calculations.  However,  at  lower  impulse 
values  (<  1  kPa  s),  the  model  overestimates  the  displacement.  This  is  likely  because  there  is 
a  threshold  force  required  for  displacement  to  occur  and  the  model  currently  does  not  take 
this  into  account. 

3.4.2  Deceleration 

For  a  qualitative  validation  of  the  deceleration  portion  of  the  model,  we  compare  model  pre¬ 
dictions  to  data  on  pedestrian  displacement  following  car  accidents.  For  forensics  purposes, 
data  has  been  collected  on  the  distance  a  pedestrian  hit  by  an  automobile  is  thrown  as  a 
function  of  the  velocity  of  the  car  at  impact  (Otte  2001).  We  assume  the  horizontal  velocity 
of  the  person  immediately  after  impact  is  equal  to  that  of  the  car  and  compare  this  data  to 
model  predictions  (Figure  3.4).  In  this  comparison,  the  model  slightly  underestimates  the 
displacement.  Because  energy  is  absorbed  upon  impact,  assuming  the  initial  velocity  of  the 
pedestrian  is  equal  to  that  of  the  car  may  not  be  accurate.  Furthermore,  the  model  may  be 
overestimating  deceleration.  The  decelerative  function  was  derived  by  scaling  animal  data 
to  humans  which  may  result  in  inaccuracies  due  to  different  tumbling  dynamics.  Also,  in 
the  data  used  for  validation,  there  are  uncertainties  in  the  car  velocity,  pedestrian  weight 


Figure  3.4:  Model  prediction  of  pedestrian  displacement  (black  line)  compared 
with  car  accident  data  (blue  circles;  Otte  2001). 
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(assumed  to  be  75  kg)  and  dynamics  of  the  impact  (e.g.,  the  height  reached  by  an  individual 
and  the  time  spent  in  contact  with  the  automobile).  Future  work  will  involve  looking  more 
critically  at  accident  data,  dummy  tests,  and  recent  derivations  on  the  relationship  between 
pedestrian  velocity  and  distance  thrown  to  improve  upon  the  translation  model  or  obtain 
better  data  for  validation  (Hague  2001). 
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4  Injury  Probit  Models 

In  this  section,  we  present  probit  models  which  relate  an  object’s  velocity  to  the  probability 
of  injury  due  to  blast.  Existing  injury  criteria  models  from  the  nuclear  effects  community  are 
reviewed  which  include  models  from  the  2013  PRCC  Report,  EM-1  Ch.  14,  and  a  presentation 
by  Dr.  John  Mercier  at  the  DTRA  HRP  meetings  in  the  early  2000s  (Department  of  the 
Army  2013;  Drake  et  al.  1993;  Reeves  2015;  Mercier  2001).  In  some  cases  where  ambiguities 
exist,  new  probit  models  are  generated. 

For  each  injury  type,  there  are  multiple  injury  levels  ranging  from  moderate  to  fatal.  In 
Section  4.1,  we  describe  the  different  injury  levels  and  contemporary  methods  for  ranking 
injury  severity.  In  Section  4.2,  the  generic  probit  model  equations  are  provided  as  well  as 
methods  for  calculating  uncertainties. 

In  Section  4.3  and  4.4,  we  review  the  existing  probit  curves  relating  the  probability  of 
tertiary  injury  to  the  maximum  velocity  obtained  by  the  human  body.  Tertiary  injury  criteria 
models  exist  for  decelerative  tumbling  and  perpendicular  impact.  Existing  probits  assume 
decelerative  tumbling  occurs  along  a  flat,  horizontal  surface  and  impact  occurs  against  a 
perpendicular,  non-deformable  surface.  Other  mechanisms  of  tertiary  injury,  such  as  impact 
against  a  deformable  surface  or  tumbling  down  a  slope,  have  not  been  thoroughly  investigated 
previously.  Based  on  the  review  of  these  probits,  we  decided  to  generate  new  probit  models 
for  injury  due  to  decelerative  tumbling  and  perpendicular  impact. 

Section  4.5  and  4.6  discuss  secondary  injuries  caused  by  penetrating  injury  and  blunt 
trauma,  respectively.  The  2013  PRCC  report  provides  injury  criteria  models  for  glass  and 
stone  penetration.  EM-1  Ch.  14  is  in  general  agreement  with  the  PRCC  report.  We  reviewed 
the  source  data  used  by  the  PRCC  report  for  probit  generation.  However,  we  were  unable 
to  re-derive  the  presented  penetration  probits.  Therefore,  using  the  same  source  data,  we 
developed  revised  penetration  probit  models.  For  blunt  trauma,  the  PRCC  report  and  EM-1 
Ch.  14  present  the  same  injury  criteria  model;  however,  we  were  unable  to  obtain  the  source 
documentation  for  this  probit. 

4.1  Injury  Severity  Levels 

The  probit  models  discussed  represent  injury  severity  levels  ranging  from  moderate  to  fatal, 
as  described  in  Table  4.1.  Traditionally,  separate  terminology  is  used  for  injury  severity 
levels  in  military  personnel  and  civilians.  For  example,  injuries  resulting  in  death  are  termed 
immediate  permanent  ineffectiveness  (IPI)  for  military  personnel  and  fatal  injury  (FI)  for 
civilians  (Department  of  the  Army  2013;  Jackson  et  al.  2013).  For  blast,  the  main  difference 


Table  4.1:  Injury  level  definitions. 


Injury  Level 

Description 

Fatal  Injury  (FI) 

Serious  Injury  (SI) 

Moderate  Injury  (MI) 
Combat  Ineffectiveness  (Cl) 

Injury  resulting  in  death 

Fractured  bones  or  ruptured  internal  organs 

Multiple  lacerations/contusions  due  to  small  missiles 
Performance  decrement  greater  than  75% 

Definitions  obtained  from  Department  of  the  Army  2013;  Jackson  et  al.  2013. 
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in  terminologies  is  that  military  personnel  are  assumed  to  be  in  a  prone  position  while 
civilians  are  assumed  to  be  standing  upon  arrival  of  the  blast  wave.  Thus,  the  velocity 
obtained  by  a  military  personnel  is  less  than  in  a  civilian  following  the  same  blast  wave; 
however,  the  velocities  that  ultimately  result  in  1P1  and  FI  are  the  same.  The  injury  levels 
listed  in  Table  4.1  are  those  relevant  to  the  probit  models  under  consideration  in  our  work: 
fatal  injury  (FI),  serious  injury  (SI),  moderate  injury  (MI),  and  combat  ineffectiveness  (Cl). 
According  to  the  PRCC  report,  1%,  2.5%  and  5%  Cl  correspond  to  negligible,  moderate, 
and  emergency  levels  of  risk  (Department  of  the  Army  2013). 

In  contemporary  data,  injury  levels  are  often  presented  using  the  Injury  Severity  Score 
(ISS).  The  ISS  is  a  scoring  system  used  to  assess  trauma  severity  on  a  whole  body  level 
(Baker  et  al.  1974).  The  ISS  is  based  upon  the  Abbreviated  Injury  Scale  (AIS).  The  AIS 
classifies  an  injury  based  on  its  location  in  the  body  and  severity.  There  are  nine  possible 
locations,  and  the  injury  severity  is  categorized  using  a  number  ranging  from  1  to  6,  where 
a  3  corresponds  to  a  serious  injury.  The  ISS  is  equal  to  the  squared  sum  of  the  highest  AIS 
scores  from  the  three  most  injured  regions  of  the  body.  We  classify  subjects  with  an  ISS 
score  of  9  or  greater  as  having  SI.  Thus,  in  some  cases,  SI  data  includes  subjects  with  fatal 
injuries. 

4.2  Probit  Model  Equation 

Probit  curves  relating  an  input  (e.g.,  velocity  and  energy)  to  injury  are  defined  by  the 
PRCC  report,  EM-1  Ch.  14  and  the  HRP  panel  using  the  following  relation  (Drake  et  al. 
1993;  Reeves  2015;  Department  of  the  Army  2013;  Mercier  2001): 

$-\pi)=P  log10(^-)  (31) 

a50 

where  pi  is  the  probability  of  an  injury  type  ranging  from  0  to  1,  x  is  the  input,  /3  is  a 
dimensionless  parameter  representing  the  slope  of  the  probit  model,  X50  is  the  input  at 
which  there  is  a  50%  chance  of  injury,  and  $  is  the  cumulative  normal  distribution  defined 
as 


$(-)  =  y=  J  e  *2  dt ■  (32) 

—  OO 

In  the  generation  of  new  injury  criteria  models,  we  will  continue  to  use  this  same  probit 
function. 

The  data  used  to  generate  probit  models  for  estimating  probability  of  injury  can  be 
represented  with  a  binomial  distribution  where  there  are  n  samples  and  the  probability  of 
injury  is  pj.  The  following  formula  defines  the  uncertainty  associated  with  the  predicted 
value  of  pp. 

2  pi{l-pi) 

a  = - 

n 

where  cr  is  the  standard  deviation.  In  this  report,  we  use  this  formula  to  give  the  error  on 
probability  estimates. 
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For  tertiary  injury,  due  to  the  range  of  probit  models  that  will  be  discussed,  the  following 
notation  will  be  used.  will  refer  to  the  velocity  at  which  there  is  an  %  percent  probability 
of  injury  j  due  to  k.  Where  i  —  0  to  100,  j  equals  the  injury  levels  (i.e. ,  FI,  Cl,  SI,  and  MI), 
and  k  is  the  injury  type  (DT  for  decelerative  tumbling  and  PI  for  perpendicular  impact). 
Similarly,  the  slopes  of  probits  will  be  represent  as 

4.3  Decelerative  Tumbling 

There  are  several  probit  models  which  calculate  the  probability  of  injury  due  to  decelerative 
tumbling  as  a  function  of  the  maximum  velocity  obtained.  Due  to  limitations  in  data,  the 
existing  probit  models  are  based  on  extrapolations  from  animal  experiments.  EM-1  Ch.  14 
provides  probit  models  for  FI  and  SI,  and  the  2013  PRCC  report  provides  a  probit  model 
for  CL  Dr.  Mercier  presented  probit  models  for  FI,  SI,  and  MI  at  the  the  HRP  meetings. 
Due  to  ambiguities  and  assumptions  underlying  existing  probits,  new  probits  were  generated 
here  using  motorcycle  accident  data  to  represent  injury  likelihood  in  urban  environments. 
Table  4.2  provides  the  probit  equations  for  all  the  decelerative  tumbling  probit  models. 


Table  4.2:  Probit  models  for  injury  due  to  decelerative  tumbling. 


Injury 

Level 

Source 

Probit  Equation 

Data  /  Assumptions 

FI 

EM-1  Ch.  14  (Reeves 
2015) 

19.55  log10R/V5o) 
I4o=45. 11  m/s 

Assumes  5%  probability  of  FI  is  equal 
to  95%  probability  of  SI.  Assumes  slope 
is  same  as  EM-1  Ch.  14  FI  perpendicu¬ 
lar  impact  probit. 

HRP  (Mercier  2001) 

19.52  log10R/V50) 
I4o=44. 39  m/s 

No  information  on  probit  generation 
methodology  or  source  data  given. 

ARA 

2.861og10R/F50) 

I4o=40.19  m/s 

Based  on  motorcycle  accident  data 
(Hurt  et  al.  1981a) 

SI 

EM-1  Ch.  14/HRP 
(Reeves  2015) 

6.42  log10R/F5o) 
1/50=20.24  m/s 

Based  on  sheep  data  (source  data  refer¬ 
ence  unknown). 

ARA 

2.40  log10R/V5o) 
F5o=9.22  m/s 

Based  on  motorcycle  accident  data 
(Hurt  et  al.  1981a) 

MI 

HRP  (Mercier  2001) 

6.42  log10(u/I40) 
140=17.51  m/s 

No  information  on  probit  generation 
methodology  or  source  data  given. 

Cl 

PRCC  (Department 
of  the  Army  2013) 

5.40  log10(w/F50) 
140=23.16  m/s 

Derived  using  goat  data  (Anderson  et 
al.  1961).  Assumes  slope  is  the  same  as 
PRCC  Cl  perpendicular  impact  probit. 
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Table  4.3:  Verification  of  probit  model  for  Cl  due  to  decelerative  tumbling. 


Reported  by  PRCC* 

Calculated 

Velocity  — >■  1%  Cl  (Negligible  Risk) 

9.1  m/s 

8.5  m/s 

Velocity  — >  2.5%  Cl  (Moderate  Risk) 

10.7  m/s 

10.1 

m/s 

Velocity  — ►  5%  Cl  (Emergency  Risk) 

11.9  m/s 

11.6 

m/s 

*Obtained  from  Department  of  the  Army  2013. 


4.3.1  PRCC  Report 

The  PRCC  report  provides  a  probit  curve  for  Cl  due  to  decelerative  tumbling  (Table  4.2). 
The  V50 ,ci,dt  is  given  in  the  PRCC  report  (23.16  m/s),  and  the  slope  Pci,dt  was  derived3. 
To  verify  the  implementation  of  this  probit  model,  the  velocities  which  result  in  1%,  2.5% 
and  5%  Cl  were  calculated  using  the  probit  model  and  compared  to  the  velocity  values  given 
in  the  PRCC  report  (Table  4.3).  These  values  are  within  ±  0.6  m/s  of  each  other. 

There  were  two  major  assumptions  in  the  Cl  probit  model  development:  (1)  injury  data 
was  extrapolated  directly  from  animals  to  humans  (i.e. ,  the  velocity  that  led  to  33%  Cl  in 
goats  was  assumed  to  be  the  same  as  the  velocity  that  led  to  33%  Cl  in  humans)  and  (2) 
because  of  the  limited  animal  data,  the  probit  slope  for  Cl  due  to  decelerative  tumbling  was 
assumed  to  be  equal  to  the  slope  of  the  perpendicular  impact  Cl  probit  (/ 3ci,dt  =  Pci, pi)- 

The  source  data  for  the  injury  criteria  model  came  from  a  blast  wave  study  (Anderson  et 
al.  1961).  In  this  study,  animals  were  placed  in  different  orientations  and  locations  within  a 
blast  tube  and  subjected  to  a  blast  wave.  Their  subsequent  maximum  velocity  was  recorded. 
The  animals  were  classified  as  injured  if  they  were  paralyzed  or  died.  Upon  review  of  the 
primary  study,  it  became  clear  that  all  the  animals  classified  as  injured  died  within  a  few 
hours  of  blast.  This  suggests  that  the  data  from  this  study  may  be  more  applicable  for 
estimating  the  likelihood  of  mortality.  The  PRCC  report  used  this  study  to  report  a  33% 
Cl  probability  for  a  velocity  of  19.2  m/s  (range:  16.5-23.8  m/s).  However,  we  were  unable 
to  derive  these  parameters  from  the  original  source  data  (Anderson  et  al.  1961). 


4.3.2  EM-1  Ch.  14 

EM-1  Ch.  14  provides  probit  curves  for  SI  and  FI  due  to  decelerative  tumbling.  The  probit 
equation  for  SI  is  not  explicitly  stated  in  EM-1  Ch.  14  but  is  rather  shown  in  Figure  60 
(Reeves  2015).  A  probit  equation  given  in  a  LFMER  report  (Fletcher  et  al.  1975)  was  found 
to  match  the  SI  probit  curve  shown  in  EM-1  Ch.  14  (see  Table  4.2  for  equation).  EM-1 
Ch.  14  and  the  Fletcher  et  al.  1975  study  both  reference  the  Middle  North  Series  reports  as 
the  source  for  the  SI  probit  model  (Richmond  et  al.  1974a;  Richmond  et  al.  1974b).  However, 

3In  the  PRCC  report,  the  velocity  for  33%  and  50%  Cl  due  to  decelerative  tumbling  are  provided 
(U33 ,CI,dt  =  19.20  m/s  and  U50 ,ci,dt  =  23.16  m/s).  Using  this  information  and  Equation  31,  we  back- 
calculate  the  probit  model  slope: 


Pci,dt 


«»-1(0.33) 

l°glo(U33,C/,Ur/U50.C/,£)T) 


(34) 


where  ^-\0.33)  =  -0.4399132 
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upon  review  of  these  published  references,  no  decelerative  tumbling  data  were  found.  Due  to 
a  lack  of  a  clear  reference  to  the  source  data  for  generation  of  this  probit,  we  cannot  verify 
the  basis  on  which  the  probit  was  developed. 

The  FI  probit  curve  presented  in  EM-1  Ch.  14  was  derived  using  assumptions  given  in 
Drake  et  al.  1978  (see  page  5-99).  These  include  (1)  a  95%  probability  of  SI  is  equivalent 
to  a  5%  probability  of  FI,  and  (2)  the  slopes  of  the  probits  for  FI  due  to  decelerative 
tumbling  and  perpendicular  impact  are  equal  (/ 3fi,dt  =  /3fi,pi  =  19.55,  see  Section  4.4.2 
for  the  perpendicular  impact  probit  slope  source).  Using  the  SI  probit  given  in  Table  4.2,  a 
maximum  velocity  of  36.6  m/s  would  result  in  95%  probability  of  SI.  Thus,  for  the  FI  probit, 
ko5 ,fi,dt  =  36.6  m/s.  In  turn,  V50 ,fi,dt  =  45.11  m/s4.  However,  as  stated  previously,  due 
to  the  uncertainty  of  the  SI  probit,  the  derived  FI  probit  is  also  questionable. 


4.3.3  HRP 

Three  probit  models  for  decelerative  tumbling  were  presented  by  Dr.  Mercier  at  the  HRP 
meetings  in  the  early  2000s  (Mercier  2001).  These  include  models  for  FI,  SI,  and  MI.  The 
probit  parameters  for  injury  due  to  decelerative  tumbling  are  given  in  Table  4.2. 

The  SI  probit  matches  that  given  in  EM-1  Ch.  14,  and  the  FI  probit  is  very  similar  to  the 
EM-1  Ch.  14  probit.  Mercier  cited  Fletcher  et  al.  1975  and  Drake  et  al.  1978  as  references  for 
the  SI  and  FI  probits.  The  MI  probit  developed  by  Mercier  has  the  same  slope  as  the  EM-1 
Ch.  14  SI  probit  but  a  slightly  lower  velocity  leading  to  50%  injury  (17.5  m/s  compared 
with  20.2  m/s).  Mercier  reports  that  the  MI  probit  is  based  on  both  sheep  data  and  human 
jumper  data;  however,  no  source  references  are  provided.  We  note  that  human  jumper  data 
would  be  more  applicable  to  impact  injuries. 


4.3.4  AMedP-8 


The  AMedP-8(C)  NATO  planning  guide  of  2009  provides  information  on  FI  due  to  the 
decelerative  tumbling  caused  by  blast  (NATO  2009).  In  this  reference,  50%  FI  for  a  prone 
personnel  in  a  random  orientation  is  related  directly  to  yield  and  static  overpressure.  NATO 
derived  this  relationship  using  a  log-linear  interpolation  and  three  data  points  obtained  from 
Drake  et  al.  1978.  Drake  et  al.  calculated  these  data  points  using  the  look-up  tables  from 
Fletcher  et  al.  1975;  however,  the  specific  velocity  which  leads  to  the  50%  FI  rate  is  not 
provided.  Using  the  Tables  A-l  and  A-3  from  Fletcher  et  al.  1975,  the  V50  velocity  used 
must  be  between  33.5  m/s  and  43.0  m/s. 

The  relationship  between  blast  parameters  and  FI  given  in  AMedP-8  was  developed  using 
the  same  methodologies  as  EM-1  Ch.  14  and  the  PRCC  report.  However,  fewer  data  points 
were  used  (i.e. ,  only  data  points  on  the  overpressure  values  at  1,  10,  100  kT  required  to 
produce  50%  FI).  AMedP-8  does  not  provide  information  on  a  probit  slope. 


4 


where  ^(O.OS)  =  -1.644854. 


_  Uo5  ,FI,DT 

U50 ,FI,DT  -  10$-i (o.05)//3 FI,DT 


(35) 
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Table  4.4:  Motorcyclist  post-crash  motion. 


Post-Crash  Motion 

Percent  of  Motorcyclists  Examined 

Stopped  near  point  of  impact 

9.1% 

Vaulted  from  motorcycle 

26.7% 

Fell  from  motorcycle 

25.9% 

Tumbled  or  rolled 

12.8% 

Slid  to  a  stop 

11.4% 

Trapped  under  motorcycle 

8.8% 

Trapped  under  other  vehicle 

2.3% 

Struck  and  dragged  by  other  vehicle 

1.2% 

Data  obtained  from  Table  8.3.2  of  Hurt  et  al.  1981a. 


4.3.5  Probit  Generation 

Because  the  decelerative  tumbling  probits  provided  in  the  PRCC  report  and  EM-1  Ch.  14 
were  based  on  very  limited  animal  data  with  large  inherent  uncertainties,  we  initiated  a 
search  for  better  data  sources.  Motorcycle  accidents,  although  highly  uncontrolled,  provide 
a  representation  of  what  might  be  observed  in  an  urban  environment  following  exposure  to 
a  blast  wave.  However,  we  would  not  expect  this  data  to  represent  injuries  in  open  field 
environments.  Thus,  we  propose  to  use  the  probits  generated  from  motorcycle  accident  data 
to  predict  injury  only  in  urban  environments. 

We  examined  a  study  on  motorcycle  accidents  done  for  the  U.S.  Department  of  Trans¬ 
portation  in  1981  (Hurt  et  al.  1981a).  This  study  contains  data  on  900  motorcycle  accidents 
and  includes  details  on  injury  severity,  collision  speeds,  helmet  use,  and  post-crash  motion. 
Data  from  motorcycle  accidents  represent  several  different  types  of  injury  due  to  contact 
with  the  motorcycle,  other  vehicles,  or  stationary  objects.  For  884  out  of  the  900  motorcy¬ 
clists  examined,  the  motorcyclist  post-crash  motion  is  given  (Table  4.4).  76.8%  of  the  data 
came  from  motorcyclists  who  vaulted  from  motorcycle,  fell  from  motorcycle,  tumbled/rollcd, 
or  slid  to  a  stop.  40%  of  the  motorcyclists  in  the  study  were  wearing  helmets.  Table  4.5 
contains  injury  versus  crash  speed  data  derived  from  the  motorcycle  accident  study.  Injury 


Table  4.5:  Proportion  of  motorcyclists  with  SI  and  FI  as  a  function  of 
crash  velocity. 


Midpoint  (range)  of 
crash  speeds  (m/s) 

Total  Motorcyclists 

Number  (%)  with  SI* 

Number  (%)  with  FI** 

2.34  (0-4.47) 

83 

24  (28.9%) 

1  (1.2%) 

6.71  (4.92-8.94) 

349 

123  (35.2%) 

7  (2.0%) 

11.18  (9.39-13.41) 

270 

166  (61.5%) 

9  (3.3%) 

15.65  (13.86-17.88) 

129 

92  (71.3%) 

18  (14.0%) 

20.12  (18.33-22.35) 

32 

26  (81.3%) 

6  (18.8%) 

24.59  (22.80-26.82) 

24 

17  (70-8%) 

8  (33.3%) 

29.06  (27.27-31.29) 

9 

7  (77.8%) 

4  (44.4%) 

33.53  (31.74-35.76) 

2 

2  (100%) 

0  (0%) 

*Derived  from  Table  8.8.6  from  Hurt  et  al.  1981a  where  SI  is  equivalent  to  an  SS>6. 
** Derived  from  Table  8.8.3  and  Table  8.8.5  from  Hurt  et  al.  1981b. 
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Figure  4.1:  Comparison  of  motorcycle  probits  for  FI  and  SI  with  data  (Hurt 
et  al.  1981a).  Goat  data  derived  from  Anderson  et  al.  1961. 


levels  were  quantified  in  the  study  using  a  “Severities  Sum”  (SS)  score,  where  an  SS  score 
of  6-12  corresponds  to  an  AIS  of  3  (Hurt  et  al.  1981a).  Thus,  we  classify  subjects  with  SS 
scores  of  6  or  greater  as  seriously  injured. 

Probits  for  FI  and  SI  as  a  function  of  velocity  were  generated.  The  probit  models  were 
parameterized  using  the  data  given  in  Table  4.5.  However,  we  excluded  the  lowest  velocity 
group  (<4.47  m/s)  because  20%  of  these  accidents  occurred  at  a  crash  velocity  of  0  m/s. 
This  indicates  that  other  injury  types  besides  decelerative  tumbling  significantly  contributed 
to  the  injuries  observed  in  this  group  (e.g.,  crushed  by  motorcycle  and  struck  by  another 
vehicle).  While  injuries  unrelated  to  decelerative  tumbling  may  be  included  in  other  velocity 
groups,  we  do  not  anticipate  that  they  impact  the  injury  curve  as  dramatically  as  expected 
in  the  slower  velocity  group.  Because  individual  crash  speeds  were  not  given,  we  used  the 
midpoint  of  the  velocity  range  for  each  group  as  the  crash  velocity.  Figure  4.1  shows  the 
probits  overlaid  on  the  accident  data  along  with  the  95%  confidence  interval.  The  data 
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Figure  4.2:  Decelerative  tumbling  probits. 


point  derived  from  the  decelerative  tumbling  study  in  goats  is  also  shown  (Anderson  et  al. 
1961).  In  this  study,  15  goats  obtained  a  velocity  between  16.6  and  20.6  m/s  with  an  average 
velocity  of  18.7  m/s.  Of  these  goats,  4  died,  leading  to  mortality  rate  of  27%5.  The  FI  probit 
generated  from  motorcyclist  data  is  consistent  with  the  data  point  from  the  goat  tumbling 
study. 

4.3.6  Analysis 

Figure  4.2  shows  a  comparison  of  all  the  decelerative  tumbling  injury  criteria  models.  The 
V50  for  FI  obtained  using  motorcyclist  data  is  between  the  V50  for  FI  from  EM-1  Ch.  14 
and  the  V50  for  Cl  from  the  PRCC  report  (see  Table  4.2).  The  PRCC  Cl  probit  model  is 

5In  the  Anderson  study,  one  of  the  goats  obtained  a  velocity  of  23.7  m/s.  However,  because  the  4  goats 
who  died  had  a  mean  velocity  of  18.0  m/s,  we  were  able  to  determine  that  the  goat  at  23.7  m/s  did  not  die 
and  exclude  it,  thereby  limiting  the  velocity  range. 
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comparable  to  the  FI  probit  models  since  all  the  goats  in  the  source  data  classified  as  injured 
died  (Anderson  et  al.  1961).  One  report  uses  the  FI  probits  from  EM-1  Ch.  14  and  the  PRCC 
report  to  estimate  a  V50  of  34.7  m/s  (Drake  et  al.  1978),  which  falls  into  the  95%  confidence 
interval  predicted  by  the  probit  derived  using  motorcycle  accident  data  (see  Figure  4.1). 
The  V50  for  SI,  obtained  using  motorcyclist  data,  is  much  less  than  that  predicted  by  EM-1 
Ch.  14  (see  Table  4.2).  The  SI  probit  from  EM-1  Ch.  14  is  based  on  sheep  data  and  Drake 
et  al.  1978  notes  that  sheep  are  a  round  animal  which  tend  to  roll  upon  impact  with  the 
surface.  This  suggests  sheep  can  withstand  a  higher  velocity  than  humans  prior  to  injury. 
Thus,  a  lower  V50  than  that  predicted  with  sheep  data  is  not  unrealistic. 

The  probit  slopes  obtained  with  the  motorcycle  accident  data  are  significantly  less  than 
the  slopes  of  the  EM-1  Ch.  14,  PRCC,  and  HRP  probit  models  (see  Figure  4.2).  When  two 
probit  models  have  different  slopes,  this  inherently  implies  that  there  is  a  threshold  velocity 
at  which  the  probits  cross.  For  example,  comparing  the  motorcycle  probits  with  the  EM-1 
Ch.  14  probits,  this  threshold  velocity  is  equal  to  32.37  m/s  for  SI  and  46.01  111/s  for  FI.  Below 
this  threshold,  the  motorcycle  probits  predict  a  higher  probability  of  injury  than  the  EM-1 
Ch.  14  probits.  Above  this  threshold,  the  motorcycle  probits  predict  a  lower  probability  of 
injury  than  the  EM-1  Ch.  14  probits.  The  EM-1  Ch.  14,  PRCC,  and  HRP  probit  models  are 
meant  to  predict  injury  due  to  decelerative  tumbling  in  open  field  environments.  Meanwhile, 
the  motorcycle  accident  data  represents  a  more  complex  environment  due  to  interactions  of 
the  body  with  the  motorcycle  and  other  objects.  Thus,  the  motorcyclist  data  are  more 
applicable  for  urban  scenarios.  Due  to  the  complex  nature  of  urban  environments,  a  higher 
probability  of  injury  at  lower  velocities  would  be  expected.  However,  at  higher  velocities, 
one  would  not  expect  the  likelihood  of  injury  to  be  less  than  that  observed  in  the  open 
field.  The  source  of  this  discrepancy  at  higher  velocities  could  arise  from  two  sources:  errors 
in  the  slopes  of  the  open  field  probits  or  the  limitations  of  using  a  probit  type  model  in 
representing  a  more  complex  injury  mechanism.  Regarding  potential  errors  in  the  slopes  of 
the  open  field  probits,  the  PRCC  Cl  and  the  EM-1  Ch.  14  FI  probit  slopes  were  assumed 
to  equal  that  of  the  perpendicular  impact  probits.  However,  during  decelerative  tumbling, 
multiple  impacts  occur,  and  the  slope  of  the  probit  represents  this  physical  process  (Drake 
et  al.  1978).  Assuming  perpendicular  and  decelerative  probit  slopes  are  equivalent  may  not 
be  an  accurate  assumption. 

To  predict  injury  due  to  decelerative  tumbling  in  open  field  environments,  we  recommend 
implementing  the  probits  from  EM-1  Ch.  14  into  HENRE  2.0.  The  EM-1  Ch.  14  probits 
were  chosen  instead  of  the  probits  presented  by  Dr.  Mercier  because  they  have  been  pub¬ 
lished  and  are,  thus,  better  documented.  All  of  the  preexisting  decelerative  tumbling  probits 
reviewed  here  have  similar  issues  regarding  the  availability  of  data,  details  on  methodologies, 
and/or  quality  of  data  used  in  their  development.  The  animals  used  to  generate  the  probit 
models  (i.e. ,  goats  and  sheep)  are  quadrupedic,  meaning  the  mechanics  of  how  they  tumble 
is  different  than  that  in  bipeds  (i.e.,  humans).  In  some  experiments,  the  animals  used  were 
anesthetized  which  likely  affects  the  tumbling  dynamics.  Finally,  the  animal  data  used  is 
very  limited  and  does  not  involve  a  wide  range  of  velocities  which  leads  to  large  uncertainties 
in  the  probit  slopes. 

To  predict  injury  in  urban  environments,  we  recommend  merging  the  motorcycle  accident 
probits  with  the  EM-1  Ch.  14  probits.  Motorcycle  accidents,  although  containing  other 
injury  types  and  in  some  cases  involving  head  protection  and  medical  treatment,  are  more 
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representative  of  the  chaotic  state  that  would  occur  following  blast  in  an  urban  environment. 
However,  motorcyclist  data  was  only  available  up  to  about  30  m/s.  At  higher  velocities,  the 
predicted  injury  from  EM-1  Ch.  14  is  greater  than  that  predicted  using  motorcycle  accident 
data.  Although  the  motorcycle  accident  probits  show  a  good  match  to  the  available  data 
(see  Figure  4.1),  it  is  possible  that  at  higher  velocities  the  probit  models  do  not  adequately 
predict  the  injury  response.  Specifically,  at  a  threshold  velocity,  the  probability  of  injury 
likely  begins  increasing  more  rapidly;  however,  a  simple  probit  model  is  not  able  to  capture 
this  effect.  Thus,  we  propose  using  the  motorcycle  accident  probit  at  low  velocities  and  the 
EM-1  Ch.  14  probit  at  high  velocities  to  predict  injury  in  urban  environments.  The  transition 
between  the  two  occurs  where  the  probits  cross.  This  threshold  velocity  is  v  =  46.01  m/s 
for  FI  and  v  =  32.37  m/s  for  SI.  Due  to  limited  data,  we  believe  the  EM-1  Ch.  14  probit 
provides  the  best  representation  of  this  transition. 

The  finalized  models  for  FI  and  SI  in  both  urban  and  open  held  environments  will  be 
implemented  in  HENRE  2.0.  In  order  to  avoid  any  discrepancies,  we  do  not  recommend  the 
inclusion  of  probits  for  other  injury  types  at  this  point  (i.e.,  MI  and  Cl). 
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4.4  Perpendicular  Impact 

In  this  section,  we  analyze  existing  probits  relating  impact  velocity  to  injury  caused  by 
perpendicular  impact  with  a  nonyielding  surface.  Again  probits  from  EM-1  Ch.  14,  the  2013 
PRCC  report,  and  PDCALC  v8.1  were  analyzed  and  are  discussed  below.  Furthermore,  due 
to  ambiguities  with  the  data,  we  again  generated  new  probit  models  using  more  contemporary 
data.  Table  4.6  provides  a  complete  list  of  probit  models  for  perpendicular  impact. 

4.4.1  PRCC  Report 

The  PRCC  report  used  data  from  a  1965  study  on  human  jumpers  to  determine  the  likelihood 
of  FI  as  a  function  of  impact  velocity  (Lewis  et  al.  1965).  This  study  provides  survival  data 
for  53  subjects  who  jumped  from  heights  ranging  from  3  to  16  stories.  The  probit  indicates 
a  V50 ,fi,pi  °f  18.08  m/s  with  a  Pfi,pi  of  5.40.  The  PRCC  report  assumes  (1)  translational 
velocities  associated  with  1%  mortality  correspond  to  50%  Cl  and  (2)  the  Cl  probit  model 
slope  is  equal  to  that  of  the  FI  probit  model  (Pci, pi  —  Pfi,pi)-  Using  this  information  the 
probit  model  for  Cl  due  to  perpendicular  impact  is  fully  specified  (see  Table  4.6). 

To  verify  the  probit  model,  the  velocities  which  result  in  1%,  2.5%,  and  5%  Cl  due  to 
impact  with  a  perpendicular,  nonyielding  surface  were  calculated  using  the  probit  model  and 
compared  to  the  velocity  values  given  in  the  PRCC  report  (Table  4.7).  These  values  are  all 

Table  4.6:  Probit  models  for  injury  due  to  perpendicular  impact. 


Injury 

Level 

Source 

Probit  Equation 

Data/  Assumptions 

FI 

EM-1  Ch.  14  (Drake 
et  al.  1993) 

19.55  log10(t;/U5o) 
V5o=  10.70  m/s 

Based  on  sheep  data  (Richmond  et 
al.  1974a).  Source  data  could  not  be 
found. 

PRCC  (Department 
of  the  Army  2013) 

5.40  log10R/F50) 
V5o=  18.08  m/s 

Based  on  human  data  (Lewis  et  al. 
1965). 

HRP  (Mercier  2001) 

13.29  log10(u/VBo) 
V50=  10.82  m/s 

No  information  on  probit  generation 
methodology  or  source  data  given. 

ARA 

7.191og10R/V50) 
V50=  15.35  m/s 

Based  on  human  data  (Lapostolle  et  al. 
2005). 

SI 

EM-1  Ch.  14  (Drake 
et  al.  1993) 

6.21  iogioR/Uso) 
1/50=  4.69  m/s 

Based  on  sheep  data  (Richmond  et  al. 
1974a). 

HRP  (Mercier  2001) 

6-211og10R/F5o) 

V50  =  5.49  m/s 

No  information  on  probit  generation 
methodology  or  source  data  given. 

ARA 

4.84  log10(i>/V5o) 
V5o=  8.61  m/s 

Based  on  human  data  (Beale  et  al.  2000; 
Dickinson  et  al.  2012). 

MI 

HRP  (Mercier  2001) 

5.931og10R/P50) 

V50  =  4.88  m/s 

No  information  on  probit  generation 
methodology  or  source  data  given. 

Cl 

PRCC  (Department 
of  the  Army  2013) 

5.40  log10(w/P5o) 
1/50=  6.71  m/s 

Assumes  1%  mortality  corresponds  to 
50%  Cl  and  the  probit  slope  is  the  same 
as  the  PRCC  FI  impact  probit  slope. 
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Table  4.7:  Verification  of  probit  model  for  Cl  due  to  perpendicular  impact. 


Reported  by  PRCC* 

Calculated 

Velocity  — >■  1%  Cl  (Negligible  Risk) 

2.6  m/s 

2.5  m/s 

Velocity  — >  2.5%  Cl  (Moderate  Risk) 

3.0  m/s 

2.9  m/s 

Velocity  — ►  5%  Cl  (Emergency  Risk) 

3.4  m/s 

3.4  m/s 

*Obtained  from  Department  of  the  Army  2013. 


within  ±0.1  m/s. 

4.4.2  EM-1  Ch.  14 

EM-1  Ch.  14  presents  probits  for  FI  and  SI  due  to  perpendicular  impact,  citing  Richmond 
et  al.  1974b  as  the  source  (Drake  et  al.  1993).  Although  this  source  did  not  include  injury 
criteria,  data  from  a  companion  document  (Richmond  et  ah  1974a)  provided  probits  for 
FI  and  SI  (see  Table  4.6).  The  probits  were  calculated  using  studies  in  which  sheep  were 
dropped  from  heights  of  0.30  to  7.6  m  on  to  a  concrete  surface.  Animals  were  dropped  in 
five  different  orientations:  head  down,  tail  down,  prone,  supine,  and  right  side  down.  SI  was 
defined  as  the  occurrence  of  fractures,  severe  organ  rupture  (liver,  heart  or  major  vessel)  and 
death.  An  impact  velocity  of  4.9  m/s  led  to  a  50%  incidence  of  SI,  and  an  impact  velocity 
of  8.11  m/s  led  to  a  1%  incidence  of  FI. 

4.4.3  HRP 

At  the  HRP  meetings,  Dr.  Mercier  presented  a  probit  model  for  FI  due  to  perpendicular 
impact  (Mercier  2001).  Similar  to  the  2013  PRCC  report,  this  probit  was  derived  with 
human  jumper  data  from  Lewis  et  ah  1965.  However,  Dr.  Mercier  made  manual  adjustments 
to  the  derived  probit  by  reducing  the  V50  by  4.57  m/s  (15  ft/s)  and  increasing  the  slope 
by  48%.  The  V50  was  reduced  to  account  for  a  random  orientation  at  impact  rather  than  a 
feet  first  orientation,  with  the  underlying  assumption  that  the  Lewis  et  ah  study  represents 
predominately  feet  first  falls.  No  justification  for  the  choice  of  4.57  m/s  was  given.  The 
increased  slope  was  based  on  sheep  data,  but,  again,  the  source  of  the  sheep  data  and  the 
reason  for  choosing  the  48%  quantity  was  not  provided.  The  finalized  probit  is  given  in 
Table  4.6. 

We  can  back-calculate  to  examine  the  probit  obtained  by  Dr.  Mercier  using  only  the 
human  jumper  data  (i.e.,  add  4.57  m/s  to  the  V50  and  divide  the  slope  by  1.48): 

<D-1(p)  =  8.9791og10(-^)  (36) 

where  V50  =  15.39  m/s.  Since  this  probit  and  the  2013  PRCC  probit  for  FI  were  derived 
with  the  same  data,  we  would  expect  them  to  match;  however,  this  is  not  the  case.  This 
discrepancy  exists  because  Dr.  Mercier  only  used  data  points  that  represented  greater  than 
10  subjects  while  the  PRCC  report  likely  used  all  data  points.  Table  4.8  gives  the  velocity 
and  mortality  data  for  jumpers  where  the  mortality  for  the  height  fallen  was  greater  than  0% 
and  less  than  100%.  Mercier  used  only  data  from  jumpers  who  fell  from  3,  4,  and  6  stories 
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Table  4.8:  Velocity  and  mortality  data  for  human  falls  from  heights. 


Height  in  Stories 

Velocity  (nr/s) 

Mortality  % 

Number 

3 

13.39 

33 

15 

4 

16.40 

64 

11 

5 

18.93 

22 

7 

6 

21.17 

91 

11 

7 

23.13 

50 

2 

Data  obtained  from  Lewis  et  al.  1965. 


(Mercier  2001).  The  inclusion  of  data  front  jumpers  from  5  and  7  stories  would  result  in  a 
reduced  probit  slope  as  seen  in  the  PRCC  probit  for  FI  due  to  impact. 

At  the  HRP  meetings,  Dr.  Mercier  also  presented  probits  relating  impact  velocity  to 
SI  and  MI.  Dr.  Mercier’s  presentation  states  these  probits  are  based  on  sheep  and  human 
jumper  data,  but  the  derivation  is  not  given.  The  SI  probit  for  impact  has  the  same  slope 
as  that  given  by  EM-1  Ch.  14;  however,  the  V50  is  greater. 

4.4.4  Probit  Generation 

Because  the  probits  reviewed  were  based  on  sheep  data  or  obtained  from  studies  on  human 
falls  with  limited  details,  we  sought  to  generate  probits  for  SI  and  FI  due  to  perpendicular 
impact  using  contemporary  free  fall  studies.  Only  data  from  studies  that  included  adults, 
out-of-hospital  mortality,  and  both  fatal  and  non-fatal  accidents  were  considered  in  the 
generation  of  revised  impact  probit  models.  Many  studies  provide  mortality  or  injury  as 
a  function  of  height.  To  determine  the  corresponding  impact  velocity,  we  use  the  formula: 
v  =  y/2 hg  where  h  is  the  height  and  g  is  the  acceleration  due  to  gravity  (9.8  m/s2). 

A  study  on  free  fall  victims  in  France  was  used  to  generate  a  probit  model  for  FI  due 
to  impact  (Lapostollc  et  al.  2005).  This  study  included  patients  who  fell  from  a  height  of 
greater  than  3  m  and  excluded  patients  less  than  12  years  of  age.  A  total  of  287  patients  were 
included  in  the  analysis.  This  study  was  chosen  to  generate  the  probit  because  it  contained 
information  on  many  aspects  of  the  fall,  including  the  circumstances  of  the  fall  (i.e.,  suicide 
attempt,  accident,  or  escape),  whether  the  impact  surface  was  hard  or  soft,  if  there  was 
preliminary  impact  before  final  impact,  and  the  part  of  the  body  that  touched  the  ground 
first  (i.e.,  head,  lower  extremities,  posterior,  anterior,  or  lateral).  For  our  purposes,  the  ideal 
dataset  would  involve  a  hard  surface,  there  would  be  no  preliminary  impact  before  final 
impact,  and  an  even  distribution  of  parts  of  the  body  would  touch  the  ground  first.  In  the 
study,  72%  of  the  impacts  were  against  a  hard  surface,  and  in  83%  of  the  cases  there  was  no 
preliminary  impact  before  the  final  impact.  Of  the  patients  included  in  the  study,  25%  were 
head  first,  24%  were  feet  first,  and  37%  had  a  side  impact  (posterior,  anterior,  or  lateral). 
In  the  remaining  subjects,  the  orientation  was  unknown.  Other  studies  did  not  provide  as 
detailed  information  on  the  circumstances  of  the  falls  (Lewis  et  al.  1965;  Beale  et  al.  2000) 
and,  therefore,  were  not  used  in  the  FI  probit  generation. 
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Figure  4.3:  Probit  models  for  FI  due  to  perpendicular  impact  compared  with 
data. 


Figure  4.3  shows  the  FI  probit  and  the  95%  confidence  intervals  along  with  the  data  from 
the  Lapostollc  et  al.  study  used  to  generate  the  probit.  For  comparison,  data  from  other 
free  fall  studies  are  shown  (Lewis  et  al.  1965;  Beale  et  al.  2000)  along  with  existing  probits. 
The  probit  model  shows  a  reasonable  fit  to  the  data  that  was  used  in  the  generation  of  the 
PRCC  probit  model  (Lewis  et  al.  1965).  The  data  from  Beale  et  al.  suggests  the  probit 
model  should  have  a  smaller  slope. 
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Figure  4.4:  Generated  SI  probit  model  for  perpendicular  impact  compared  with 
data. 


For  generating  a  probit  for  SI  due  to  impact,  we  again  examined  free  fall  studies.  Al¬ 
though  the  study  used  for  the  FI  probit  generation  contained  more  detailed  information,  it 
did  not  include  injury  severity  data  as  a  function  of  height.  Two  other  studies,  comparable 
in  quality,  were  found  which  supplied  individual  ISS  data  for  both  fatal  and  non-fatal  in¬ 
juries,  including  patients  who  died  out-of-hospital  (Beale  et  al.  2000;  Dickinson  et  al.  2012). 
Figure  4.4  shows  the  probit  model  obtained  through  analyzing  the  individual  data. 

4.4.5  Analysis 

Figure  4.5  compares  all  the  FI  perpendicular  impact  injury  criteria  models,  and  Figure  4.6 
compares  all  the  SI,  MI,  and  Cl  perpendicular  impact  injury  criteria  models.  The  probits 
presented  by  EM-1  Ch.  14  are  based  on  sheep  data,  while  the  PRCC  probit  and  the  probits 
generated  here  are  based  on  human  data.  Due  to  the  large  physical  differences  between 
human  and  sheep,  it  is  likely  that  the  same  impact  velocity  would  result  in  different  injuries. 
Thus,  it  is  best  to  use  human  data,  if  available,  for  generating  the  probit.  However,  one 
criticism  of  the  human  data  used  is  that  the  orientation  of  jumpers  would  likely  result  in  a 
feet  first  landing  (Drake  et  al.  1978),  and  following  a  nuclear  blast  scenario  the  orientation 
of  body  during  impact  would  be  more  random.  The  study  used  by  the  PRCC  report  does 
not  provide  information  on  the  orientation  of  the  subjects  upon  impact  (Lewis  et  al.  1965). 
If  impact  was  predominately  feet  first,  the  reported  V50  value  for  FI  of  18.08  m/s  is  likely  too 
high  because  the  legs  would  allow  the  body  to  decelerate  before  major  organs  or  the  head 
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Figure  4.5:  Probit  models  for  FI  due  to  perpendicular  impact. 


impacted  the  ground. 

Rather  than  revert  to  sheep  data,  we  decided  to  search  the  literature  for  more  recent 
studies  on  human  falls  to  generate  revised  probit  models.  The  probit  generated  here  for  FI 
due  to  impact  falls  between  the  probits  given  by  EM-1  Ch.  14  and  the  PRCC  report  (see 
Figure  4.5).  The  V50  is  greater  than  that  predicted  by  EM-1  Ch.  14  and  less  than  that 
predicted  by  the  PRCC  report  (Table  4.6).  The  slope  (3  is  greater  than  that  from  the  PRCC 
report  and  less  than  that  from  EM-1  Ch.  14. 
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Impact  Velocity  (m/s) 

Figure  4.6:  Probit  models  for  SI,  Cl  and  MI  due  to  perpendicular  impact. 


A  comparison  of  the  derived  SI  probit  with  existing  SI  and  Cl  probit  models  is  shown 
in  Figure  4.6.  The  probit  derived  here  has  a  larger  V50  value  than  the  other  probit  models, 
but  a  very  similar  slope.  This  results  in  greater  SI  predictions.  The  Cl  probit  from  PRCC  is 
very  similar  to  the  derived  probit,  while  the  SI  probit  from  EM-1  Ch.  14,  which  was  derived 
using  sheep  data,  is  the  least  similar.  Of  the  reviewed  SI  probits,  the  generated  SI  probit 
model  is  the  only  probit  based  entirely  on  human  data.  However,  this  probit  was  generated 
using  digitized  data.  It  is  possible  that  there  were  overlapping  data  points  which  could  not 
be  distinguished  and,  therefore,  were  classified  as  only  one  subject. 

We  recommend  the  inclusion  of  the  SI  and  FI  probits  generated  here  into  HENRE  2.0. 
The  FI  probit  was  generated  using  data  that  provided  detail  on  subject  landing  orientation. 
For  SI,  Cl,  and  MI,  the  probit  generated  here  is  the  only  one  based  entirely  on  human  data. 
Although  probits  for  Cl  and  MI  exist,  we  recommended  not  including  those  injury  levels.  The 
Cl  probit  was  generated  by  assuming  1%  mortality  corresponded  to  50%  Cl  (Department 
of  the  Army  2013),  and  no  source  documentation  or  description  of  the  MI  probit  generation 
process  is  available  (Mercier  2001). 
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4.5  Penetrating  Debris 

Both  EM-1  Ch.  14  and  the  PRCC  report  discuss  blast  injuries  caused  by  penetrating  debris. 
Developing  injury  criteria  for  penetrating  debris  is  challenging  because  the  likelihood  of 
injury  is  dependent  on  several  factors  including  the  missile  impacting  area,  shape,  density, 
and  velocity.  Furthermore,  in  order  to  relate  blast  to  penetrating  injury,  it  is  necessary  to 
predict  the  particle  shape  and  velocity  distribution  that  an  individual  will  be  subjected  to. 
A  summary  of  the  existing  probit  models  and  models  generated  here  is  given  in  Table  4.9. 
Each  of  these  probits  will  be  discussed  here. 

4.5.1  EM-1  Ch.  14/PRCC 

Here,  we  are  first  interested  in  examining  existing  probit  models  in  EM-1  Ch.  14  and  the 
PRCC  report  which  relate  the  likelihood  of  penetration  to  missile  characteristics.  Both 
documents  discuss  penetration  injury  caused  by  glass  and  non-glass  objects.  EM-1  Ch.  14 

Table  4.9:  Probit  models  for  injury  due  to  missile  penetration. 


Injury 

Level 

Source 

Probit  Equation 

Data/  Assumptions 

SI  (glass) 

ARA 

2.851og10(^) 

F50  =  89.45 

Based  on  canine  data  (Bowen  et  al. 
1956).  Assumes  the  likelihood  of  SI 
equals  the  likelihood  that  a  glass  shard 
will  penetrate  the  abdominal  wall. 

MI  (glass) 

ARA 

2.851og10(^) 

F50  =  44.93^ChC 

Based  on  canine  data  (Bowen  et  al. 
1956).  Assumes  the  likelihood  of  MI 
equals  the  likelihood  that  a  glass  shard 
will  penetrate  the  skin. 

MI  (stone) 

ARA 

2.851og10(^) 

F50  =  178.5^^ 

Based  on  goat/human  data  (Sperrazza 
et  al.  1967).  Assumes  the  probit  slope 
is  equal  to  the  probit  slope  of  the  glass 
penetration  probit  and  the  likelihood  of 
MI  equals  the  likelihood  that  a  glass 
shard  will  penetrate  the  skin. 

Cl  (glass) 

PRCC  (Department 
of  the  Army  2013) 

9'1  log10(yP-) 

Ho  =  128  f  for  O.lg 
Ho  =  75  =  for  l.Og 
Ho  =  58“  for  lOg 

Based  on  canine  data  (Bowen  et  al. 
1956).  Assumes  the  likelihood  of  Cl 
equals  the  likelihood  a  glass  shard  will 
penetrate  the  abdomen.  We  were  un¬ 
able  to  reproduce  this  probit  using  the 
source  data. 

Cl  (stone) 

PRCC  (Department 
of  the  Army  2013) 

9-Hog10(*) 

Ho  =  188  f  for  O.lg 
Ho  =  104=  for  l.Og 
Ho  =  61=  for  lOg 

Based  on  goat/human  data  (Sperrazza 
et  al.  1967).  Assumes  the  probit  slope 
is  equal  to  the  probit  slope  of  the  glass 
penetration  probit  and  the  likelihood  of 
Cl  equals  the  likelihood  stone  will  pene¬ 
trate  the  skin.  We  were  unable  to  repro¬ 
duce  this  probit  using  the  source  data. 
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Table  4.10:  The  50%  penetration  velocities  of  steel  and  stone. 


Missile 

Stone  Mass  (grams) 

Ballistic  limits 

Bare  Skin 

Military  Uniform 

Steel  sphere* 

1.0 

60  m/s  (197  ft/s) 

154  m/s  (505  ft/s) 

2.0 

52  m/s  (171  ft/s) 

137  m/s  (450  ft/s) 

10 

40  m/s  (131  ft/s) 

111  m/s  (364  ft/s) 

Stone** 

0.1 

188  m/s  (617  ft/s) 

408  m/s  (1340  ft/s) 

1.0 

104  m/s  (340  ft/s) 

354  m/s  (1160  ft/s) 

10 

61.0  m/s  (200  ft/s) 

171  m/s  (560  ft/s) 

*Values  obtained  from  Sperrazza  et  al.  1967. 

**Values  obtained  from  Department  of  the  Army  2013. 


discusses  a  British  study  which  examined  the  effect  of  fragment  shape  on  penetration  velocity 
for  a  53  mg  missile  (UK  Ministry  of  Defense  1969) 6 .  It  was  found  that  particles  shaped  as 
spheres  required  a  lower  velocity  than  randomly  shaped  fragments  for  penetration.  Although 
providing  useful  qualitative  information,  the  missile  sizes  and  velocities  explored  in  this  study 
were  deemed  irrelevant  for  nuclear  effects  data  (Reeves  2015). 

A  second  study,  discussed  by  both  EM-1  Ch.  14  and  the  PRCC  report,  calculated  the 
ballistic  limits  of  isolated  goat  and  human  skin  exposed  to  steel  fragments  weighing  up  to 
15  g  (Sperrazza  et  ah  1967).  This  study  provides  the  following  relationship  between  the 
ballistic  limit  (V50)  and  the  fragment  mass 

U5o  =  K-  +  b  (37) 

m 

where  for  skin  K  =  1.2471  •  106  g/(m  s),  b  =  22.03  m/s,  a  is  the  cross  sectional  area  of  the 
object,  and  m  is  the  mass.  For  military  uniforms,  K  =  2.6107  •  106  g/(m  s)  and  b  =  73.51 
m/s.  The  ballistic  limits  for  1.0,  2.0,  and  10  g  steel  spheres  are  given  in  Table  4.10.  These 
values  correspond  to  a  steel  density  of  8  gm/cm3.  In  order  to  calculate  the  ballistic  limits  for 
stone,  the  PRCC  report  uses  Equation  37  and  assumes  the  drag  coefficient  of  the  stones  is 
1.0.  Thus,  the  value  of  a/m  is  equal  to  the  acceleration  coefficients  of  the  stones  which  were 
previously  determined  (Department  of  the  Army  2013).  Table  4.10  provides  the  resulting 
ballistic  limits  for  stones7. 

EM-1  Ch.  14  does  not  provide  a  probit  slope  for  stone  or  steel  penetration,  while  the 
PRCC  report  assumes  the  stone  penetration  probit  slope  is  equivalent  to  the  glass  penetration 
probit  slope.  Based  on  the  digitization  of  Figure  D-24  from  the  PRCC  report  this  slope  is 
approximately  9.4. 

For  glass  penetration,  the  PRCC  report  and  EM-1  Ch.  14  reference  a  study  on  glass 
abdominal  penetration  in  dogs  (Bowen  et  ah  1956).  EM-1  Ch.  14  reports  on  another,  more 
complex,  study  which  examined  the  effects  of  penetration  angle,  fragment  mass  and  impact 
velocity  on  V50  (Fletcher  et  ah  1980).  In  this  study,  experiments  were  performed  in  which 

6EM-1  Ch.  14  states  the  missiles  were  58  mg,  while  the  UK  source  document  states  they  were  53  mg. 

'These  values  were  pulled  directly  from  the  PRCC  2013  report;  however,  we  were  unable  to  reproduce 
these  values  precisely  given  the  acceleration  coefficients  provided  by  the  PRCC  report. 
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Table  4.11:  Glass  impact  velocities  required  for  abdominal  penetration  in  dogs. 


Mass  (g) 

Velocity  (m/s) 
1%  Penetration 

50%  Penetration 

99%  Penetration 

0.1 

75.0 

134 

232 

1.0 

44.0 

76.0 

131 

10 

30.7 

53.0 

91.5 

Data  derived  from  Bowen  et  al.  1956. 


animals  were  placed  in  front  of  windows  and  exposed  to  overpressure  waves.  However,  no 
probit  curves  were  derived  from  this  study. 

The  PRCC  report  provides  a  probit  curve  based  on  the  canine  abdominal  penetration 
study  (Bowen  et  al.  1956).  In  the  penetration  study,  the  threshold  velocity  for  abdominal 
penetration  is  given  as 

,  .  100.28  .  . 

log10(m)  =  —2.3054  H -  (38) 

H) 


or 


100.28 

vo 


(39) 


where  m  is  mass  of  the  missile  in  grams  and  Vq  is  the  threshold  velocity  in  m/s.  For  this 
equation  to  be  valid,  m  must  be  greater  than  mo  (mo  =  0.00495  g).  The  relationship  between 
missile  velocity  v  and  probability  of  injury  P  is  given  as 


l°g10(— )  =  0.4842P. 
v0 


(40) 


Table  4.11  provides  the  calculated  velocities  required  to  result  in  1%,  50%  and  99%  abdom- 
enal  penetration.  This  table  does  not  precisely  agree  with  the  Table  D-15  from  the  PRCC 
report  or  Table  14A-15  from  EM-1  Ch.  14,  even  though  both  were  derived  using  the  same 
study  (Bowen  et  al.  1956).  The  PRCC  report  presents  a  probit  model  determined  with  data 
from  this  table.  However,  since  we  were  unable  to  reproduce  the  table,  we  were  also  unable 
to  reproduce  the  probit  parameters. 


4.5.2  Probit  Generation 

In  this  section,  we  derive  new  penetration  injury  probits  for  stone  and  glass  because  (1)  we 
were  unable  to  reproduce  the  glass  penetration  probit  developed  in  the  PRCC  report  and  (2) 
we  wanted  to  generate  a  probit  for  a  continuous  range  of  object  masses.  The  PRCC  report 
generated  three  separate  probits  for  0.1  g,  1.0  g,  and  10  g.  Instead,  our  aim  is  to  include 
the  mass  in  the  probit  equation.  It  was  determined  that  the  probit  model  for  penetration 
should  be  of  the  following  form: 

t-1W  =  /?log1o(^)  (41) 

The  function  given  by  Equation  41  was  chosen  because  penetration  is  dependent  on  the 
kinetic  energy  of  the  missile  (\mv2).  However,  the  penetration  surface  area  likely  increases 
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with  the  mass,  causing  the  needed  energy  for  penetration  to  increase.  According  to  Equa¬ 
tion  41,  the  needed  impact  energy  for  50%  penetration  is  F^im}~a.  Thus,  the  value  for  a 
should  be  less  than  1  (i.e.,  as  mass  increases,  the  energy  needed  for  penetration  increases). 

To  determine  the  value  for  a  in  the  stone  penetration  model,  we  used  the  previously 
reported  ballistic  limit  data  for  stone  penetration  of  skin  (Table  4.10).  Since  the  stone  data 
is  for  50%  penetration,  Equation  41  can  be  reduced  to 

F50  =  \mav2  (42) 


Least  squares  regression  of  Equation  42  to  the  data  in  Table  4.10  resulted  in  an  a  value  equal 
to  0.5  and  an  F50  value  of  178.5  m2/s2  kg1/2.  Thus,  the  following  equation  defines  the  stone 
penetration  probit  model 


<f>_1 


P  logm 


/  ^  m0.5v2\ 

V  2  A50  / 


(43) 


where  /3  will  be  determined  by  assuming  the  same  slope  as  the  glass  penetration  probit. 

To  develop  probits  for  glass  penetration,  data  on  skin  and  abdominal  penetration  in  dogs 
were  used  (Bowen  et  al.  1956).  To  avoid  the  possibility  that  the  abdominal  penetration 
probability  is  ever  greater  than  the  skin  penetration  probability,  the  probit  slopes  for  the 
skin  and  abdominal  penetration  models  were  set  equal  to  each  other,  and  the  models  were 
optimized  using  least  squares  regression,  a  was  originally  included  in  the  model  optimiza¬ 
tion,  and  it  was  found  that  the  value  was  approximately  equal  to  0.5.  To  keep  the  probits 
consistent,  the  value  of  a  was  set  to  0.5  and  the  model  was  re-optimized.  The  hnal  glass 
penetration  probits  compared  with  data  are  shown  in  Figure  4.7. 


Figure  4.7:  Probit  model  for  glass  penetration  of  skin  and  abdomen.  Data 
obtained  from  a  study  in  dogs  (Bowen  et  al.  1956). 
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In  finalizing  the  stone  penetration  probit,  we  assumed  the  slope  was  equal  to  that  of  the 
glass  penetration  probit,  the  same  assumption  that  was  used  in  the  PRCC  report.  Thus,  the 
stone  penetration  probit  is: 

/  ^  m°-5v2  \ 

“lo^(2wj  (44) 

where  F50  =  178.5  m2/s2  kg1/2. 

These  probits  provide  information  on  the  likelihood  of  penetration  but  do  not  relate 
that  to  injury.  The  PRCC  report  assumes  50%  penetration  corresponds  to  50%  Cl;  however, 
there  are  more  mechanistic  ways  to  relate  penetration  likelihood  to  injury.  This  would  involve 
incorporating  environmental  information  on  the  mass  and  number  of  missiles  a  person  would 
encounter,  and  determining  the  probability  of  penetration  that  leads  to  a  specified  injury 
level.  This  is  a  needed  step  which  will  be  discussed  more  in  Section  6.  For  now,  we  use 
the  same  assumption  as  the  PRCC  report  and  equate  the  probability  of  skin  and  abdominal 
penetration  to  the  probability  of  MI  and  SI,  respectively. 

4.5.3  Analysis 

A  comparison  of  the  new  probits  with  existing  probits  is  shown  in  Figure  4.8.  For  this 


Figure  4.8:  Comparing  penetration  probits  for  1.0  gram  missiles. 
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comparison,  the  mass  of  the  missile  is  set  to  1.0  g.  To  compare  the  probits  generated  here 
to  preexisting  ones,  we  rearranged  Equation  43  to  obtain: 

J>-1  =  2/31og10f^  (45) 

where 


C50  = 


(46) 


Because  our  models  were  derived  with  the  same  data  as  the  PRCC  probit  models,  we  would 
expect  them  to  be  quite  similar.  This  is  indeed  the  case.  The  V50  from  the  MI  stone 
penetration  model  generated  here  is  equivalent  to  the  V/0  from  the  Cl  stone  penetration 
model  generated  in  the  PRCC  report.  Similarly  the  V5a  from  the  SI  glass  penetration  model 
is  equivalent  to  the  V50  from  the  Cl  glass  penetration  model  generated  by  PRCC.  The  main 
difference  lies  in  the  injury  severity  classihcation.  We  recommend  the  inclusion  of  the  re¬ 
derived  probits  into  HENRE  2.0. 


4.6  Blunt  Trauma 

Blunt  trauma  occurs  when  a  person  is  struck  by  an  object  which  has  sufficient  kinetic  energy 
to  cause  injury  without  penetration.  Like  other  injuries,  the  effects  of  blunt  trauma  are 
highly  dependent  on  the  region  of  the  body  that  is  struck;  the  head,  abdomen,  and  thorax 
are  the  most  sensitive. 

The  PRCC  report  presents  a  probit  relating  the  impact  energy  of  an  object  to  the  like¬ 
lihood  of  blunt  trauma  induced  CL  This  probit  is  based  on  data  from  a  study  done  by  the 
British  during  the  Second  World  War  (Zuckerman  et  al.  1944)8.  In  this  study,  lumps  of  hard 
crater  debris  and  softer  lumps  of  clay  were  impacted  against  small  animals  and  dried  human 
skulls.  Injury  was  recorded  if  death,  unconsciousness,  major  bone  fracture,  or  severe  rupture 
or  hemorrhage  occurred.  Thus,  although  the  PRCC  report  used  this  data  to  define  Cl,  it  is 
representative  of  SI  as  well. 

The  probit  provided  in  Figure  D-26  of  the  PRCC  report  was  digitized  to  obtain  the 
following  equation: 

$_1(p)  =  2.51  log  f 

V-^50 

where  E  is  the  energy  of  the  missile  and  E50  =  67.37  J.  EM-1  Ch.  14  uses  the  probit  model 
developed  in  the  PRCC  report  to  provide  V50  values  of  17.4,  7.6,  and  5.5  m/s  for  missile 
weights  of  0.45,  2.3,  and  4.5  kg  which  agrees  with  the  E50  value  derived  through  digitization. 
We  did  not  seek  to  generate  a  new  probit  for  blunt  trauma.  We  recommend  including  the 
probit  given  by  Equation  47  into  HENRE  2.01,  with  the  aim  of  eventually  replacing  it  with 
a  more  sophisticated  model. 


8We  were  unable  to  obtain  this  original  source.  The  information  provided  here  was  given  by  Drake  et  al. 
1978. 
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5  Summary 

This  work  outlines  a  translation  model  linked  with  velocity  probit  models  that  predict  the 
probability  of  injury  following  complex  blast  wave  exposures.  These  models  can  be  linked 
with  nuclear  environment  tools  that  predict  blast  parameters  in  urban  areas  to  estimate  the 
potential  fatalities  and  casualties  due  to  blast. 

5.1  Translation  Models 

The  translation  model  requires  the  following  inputs:  the  dynamic  pressure  wave,  the  time- 
course  of  the  wind  velocity,  and  the  acceleration  coefficient  of  the  object  of  interest.  The 
acceleration  coefficient  is  a  function  of  the  object’s  mass,  area  presented  to  the  wind,  and  drag 
coefficient.  With  these  inputs,  the  model  calculates  the  translational  profile  (i.e.,  distance, 
velocity,  and  acceleration)  as  a  function  of  time.  This  is  an  update  to  existing  methods 
because  complex  blast  waves  can  be  directly  inputted  to  determine  an  object’s  velocity. 

5.2  Finalized  Injury  Criteria  Models  for  Inclusion  in  HENRE  2.0 

Table  5.1  provides  a  summary  of  all  probit  models  recommended  for  inclusion  into  HENRE 
2.0.  For  tertiary  injury,  these  include  probits  for  FI  and  SI,  and  for  secondary  injury,  these 
include  probits  for  SI  and  MI.  These  recommendations  are  based  on  a  thorough  review  of 
the  existing  injury  criteria  models  used  by  the  nuclear  effects  community. 

For  both  FI  and  SI  due  to  decelerative  tumbling,  multiple  injury  criteria  models  exist. 
These  include  previously  generated  probits  presented  in  EM-1  Ch  14  and  by  Dr.  Mercier  at 

Table  5.1:  Summary  of  probit  models  recommended  for  inclusion  in  HENRE  2.0 


Fatal  Injury 

Serious  Injury 

Moderate  Injury 

Decelerative  Tumbling 

if  v  <  46.01  m/s: 

if  v  <  32.37  m/s: 

None 

(Urban  Environment) 

2.861og10(^); 

U50  =  40.19™ 
if  v  >  46.01  m/s: 
19.55  log10(p^); 
^50  =  45.11  = 

2.401og10(^); 

U50  =  9.22™ 
if  v  >  32.37  m/s: 
6-42  logic  (^); 

U50  =  20.24™ 

Decelerative  Tumbling 

19-55  logndpD); 

6.421og10(^);U5o  = 

None 

(Open  Field) 

U5o=  45.11  = 

20.24— 

S 

Perpendicular  Impact 

7.191og10(^); 

U50  =  15.35  = 

4.841og10(^);  U50  = 
8.61  = 

None 

Glass  Penetration 

None 

2.851og10(^); 

P50  =  89.45 

2.851og10(^); 

F50  =  44.93 

Stone  Penetration 

None 

None 

2.851og10(^); 

p.  0.5  2 

P50  =  178.5^^^^ 

Blunt  Trauma 

None 

2-511°Sio(^); 

E50  =  67.37J 

None 
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the  HRP  (Reeves  2015;  Mercier  2001)  and  probits  developed  here  using  motorcycle  accident 
data  (Hurt  et  al.  1981a).  Previously  generated  models  are  based  on  unobtainable  data 
and/or  the  model  generation  process  was  unclear.  The  motorcycle  accident  data  provides 
a  more  verihable  and  accurate  representation  of  injury  probabilities  that  might  be  observed 
in  an  urban  environment.  Thus,  we  recommend  the  inclusion  of  the  probits  derived  using 
motorcycle  data  into  HENRE  2.0  to  predict  injury  in  urban  environments  at  low  velocities. 
For  injury  predictions  in  open  held  environments,  the  probits  given  by  EM-1  Ch.  14  and 
presented  by  Dr.  Mercier  are  very  similar.  However,  more  description  was  given  by  EM-1 
Ch.  14  on  probit  generation  and,  thus,  the  probits  from  EM-1  Ch.  14  are  recommended  for 
inclusion  into  HENRE  2.0  to  predict  decelerative  tumbling  injuries  in  open  held  environments 
as  well  as  injuries  in  urban  environments  at  high  velocities  (Reeves  2015). 

For  Cl  and  MI  due  to  decelerative  tumbling,  we  did  not  generate  new  probits.  The 
existing  Cl  model  for  decelerative  tumbling  is  based  on  one  data  point  from  a  goat  study 
in  which  all  of  the  injured  goats  died  (Anderson  et  al.  1961).  Furthermore,  the  slope  was 
assumed  to  equal  that  of  the  perpendicular  impact  model  developed  by  PRCC.  The  source 
data  and  derivation  of  the  MI  model  is  not  described  and  warrants  further  investigation.  For 
these  reasons,  we  do  not  recommend  the  inclusion  of  a  Cl  or  MI  model  at  this  time. 

For  both  FI  and  SI  due  to  perpendicular  impact  with  a  non-yielding  surface,  multiple 
probits  exist.  These  include  preexisting  probits  (Department  of  the  Army  2013;  Drake  et 
al.  1993;  Mercier  2001)  and  probits  generated  here  using  data  from  human  free  fall  studies 
(Lapostolle  et  al.  2005;  Beale  et  al.  2000;  Dickinson  et  al.  2012).  These  probits  were  derived 
using  sheep  data,  human  data,  or  a  combination  of  the  two.  Due  to  assumptions  required 
to  translate  from  sheep  to  humans,  we  recommend  including  a  model  that  is  based  solely 
on  human  data.  For  FI,  these  include  the  models  derived  in  PRCC  report  and  by  us  using 
human  free  fall  studies.  Because  the  FI  probit  generated  here  used  a  study  that  contained 
more  detailed  information  on  falling  dynamics,  we  recommend  the  inclusion  of  this  probit 
into  HENRE  2.0.  For  SI,  the  probit  generated  by  us  is  the  only  one  based  on  human  data, 
and  therefore,  we  recommend  the  inclusion  of  this  probit. 

Probits  also  exist  for  Cl  and  MI  due  to  perpendicular  impact.  The  MI  probit  is  poorly 
documented,  with  no  the  source  data  description  (Mercier  2001).  The  Cl  probit  was  derived 
by  assuming  1%  mortality  corresponds  to  50%  Cl  (Department  of  the  Army  2013).  Due  to 
the  underlying  assumptions  and  lack  of  documentation,  we  do  not  recommend  inclusion  of 
the  Cl  or  MI  probits. 

Secondary  injury  criteria  models  are  available  for  Cl,  SI,  and  MI.  The  only  preexisting 
penetration  injury  criteria  models  were  found  in  the  PRCC  report;  however,  we  were  unable 
to  re-derive  the  stone  or  glass  penetration  probit  using  the  source  data.  Here,  we  developed 
new  probits,  using  the  same  source  data  as  the  PRCC  report,  that  include  the  mass  of 
the  object  in  the  probit  equation.  We  assumed  that  the  likelihood  of  a  missile  penetrating 
the  skin  was  equivalent  to  MI  and  the  likelihood  of  a  missile  penetrating  the  abdomen  was 
equivalent  to  SI.  Because  the  source  data  used  by  the  PRCC  report  did  not  contain  data  on 
the  likelihood  of  stone  penetrating  the  abdomen,  no  SI  model  for  stone  exists. 

Finally,  for  blunt  trauma,  only  a  probit  model  for  SI  exists.  We  were  unable  to  obtain 
the  source  documentation  for  this  blunt  trauma  model  (Zuckerman  et  al.  1944).  This  probit 
is  based  on  data  obtained  by  performing  impact  studies  with  small  animals  and  dried  hu¬ 
man  skulls  (Drake  et  al.  1978).  We  propose  the  inclusion  of  this  blunt  impact  model  into 
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Figure  5.1:  Final  tertiary  injury  probit  models. 


HENRE  2.0  so  that  a  complete  picture  of  injury  types  is  implemented.  In  the  future,  more 
sophisticated  models  of  blunt  injury  should  be  explored. 

Figure  5.1  compares  the  tertiary  probit  models  recommended  for  implementation  in 
HENRE  2.0.  The  models  for  SI  and  FI  predict  that,  at  low  velocities,  injuries  due  to  decel- 
erative  tumbling  in  urban  environments  are  more  likely  than  injuries  due  to  perpendicular 
impact.  In  decelerative  tumbling,  multiple  impacts  occur,  which  could  lead  to  an  increased 
likelihood  of  a  head  injury  when  compared  with  perpendicular  impact  (Drake  et  al.  1978). 
Also,  although  the  generated  models  improve  upon  the  existing  probits,  further  model  de¬ 
velopment  may  be  warranted.  Work  could  be  done  to  improve  the  probit  model  structure 
and/or  under  lying  data.  In  addition  to  impact  velocity,  other  metrics,  such  as  age  and  ori¬ 
entation,  are  used  elsewhere  to  predict  mortality  (Dickinson  et  al.  2012;  Lapostollc  et  al. 
2005).  These  models  could  be  explored  to  see  if  they  provide  a  more  accurate  representation 
of  the  data. 
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Figure  5.2:  Final  penetration  injury  probit  models  for  1.0  gram  missiles. 


The  penetrating  injury  probits  for  stone  and  glass  with  masses  of  1.0  g  are  shown  in 
Figure  5.2.  Given  the  same  missile  mass  and  velocity,  the  likelihood  of  MI  and  SI  from 
glass  penetration  is  greater  than  the  likelihood  of  SI  from  stone  penetration.  These  models 
assume  that  the  probability  of  SI  and  MI  are  equivalent  to  the  probability  of  abdominal  and 
skin  penetration,  respectively.  However,  it  is  possible  that  a  stone  penetration  of  the  skin  is 
actually  more  severe  than  a  glass  shard  penetration  due  to  a  larger  impact  area. 
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6  Conclusions  and  Future  Work 


This  work  provides  an  incremental  improvement  to  blast  injury  calculations;  however,  there 
are  several  facets  of  the  modeling  techniques  discussed  that  can  be  improved  upon  in  the 
future. 

Future  work  to  improve  the  translation  model  should  involve  refining  the  deceleration 
portion  of  the  model  and  incorporating  directional  components  for  the  wind  velocity  and 
dynamic  pressure.  Contemporary  work  has  been  done  to  estimate  the  distance  traveled  by 
a  pedestrian  who  has  been  hit  by  a  vehicle.  Models  have  been  developed  which  predict  the 
distance  a  pedestrian  is  thrown  based  on  the  vault  angle  and  impact  speed  (Evans  et  al.  1999). 
These  models  can  be  used  to  assess  or  replace  the  deceleration  portion  of  the  translation 
model  reviewed  here.  As  for  adding  directional  components,  the  translation  model  currently 
assumes  the  direction  of  the  dynamic  wave  is  constant,  but  in  urban  environments  this 
is  likely  not  the  case.  Thus,  future  work  should  involve  representing  the  blast  wave  and 
translational  profile  using  vector  notation. 

Probit  models  for  FI  and  SI  due  to  decelerative  tumbling  and  perpendicular  impact  have 
been  generated  using  data  obtained  from  human  accidents.  To  obtain  more  precise  probits, 
a  more  thorough  epidemiological  analysis  needs  to  be  performed.  For  example,  in  generating 
the  decelerative  tumbling  probit,  we  used  data  from  a  study  on  motorcycle  accidents.  This 
study  grouped  together  subjects  who  did  or  did  not  have  helmets  and  had  different  types 
of  motion  following  the  crash.  An  ideal  analysis  would  involve  using  the  raw  data  from  this 
study  to  select  subjects  with  specific  criteria  (e.g.,  only  include  subjects  whose  injuries  were 
due  to  tumbling).  A  similar  analysis  could  be  done  with  data  on  free  fall  victims.  This  would 
result  in  probits  which  are  based  on  the  most  appropriate  data  for  representing  the  injuries 
observed  following  blast  waves  in  urban  environments. 

Probit  models  were  developed  for  SI  and  MI  due  to  penetrating  injuries.  These  models 
provide  the  likelihood  that  a  given  fragment  will  penetrate  the  skin  or  abdomen.  However, 
relating  fragment  penetration  to  injury  has  not  yet  been  done.  Future  work  will  involve 
determining  this  relationship  either  by  using  past  work  or  by  developing  new  models.  The 
Operational  Requirement-based  Casualty  Assessment  (ORCA)  software  package  contains 
models  which  track  penetrating  fragments  and  determine  the  likelihood  of  injury  caused  by 
the  fragment  based  on  the  region  of  the  body  penetrated.  In  turn,  the  likelihood  of  injury 
due  to  multiple  fragments  can  be  ascertained. 

For  glass  penetration,  more  detailed  models  which  relate  the  shattering  of  a  window  to 
injury  probability  have  been  generated.  These  include  a  Shard  Fly-Out  Model  (SFOM)  and 
the  Multi-Hit  Glass  Penetration  (MHGP)  Model  (Meyer  et  al.  2004).  The  SFOM  determines 
if  and  how  a  window  breaks  following  a  blast  wave  and  generates  a  distribution  of  glass  shards 
with  details  on  the  shape,  mass,  velocity,  and  orientation  of  each  shard.  This  model  includes 
details  on  the  window  type  and  thickness.  The  MHGP  Model,  in  turn,  determines  whether 
injury  occurs  based  on  the  glass  shards  which  come  into  contact  with  an  individual.  The 
injury  criteria  is  based  on  the  depth  of  penetration  of  glass  shards  which  is  given  by 

777/ 

Xmax  =  ln(l  +  0.001216E)  (48) 

O.loOoAp 

where  xmax  is  the  maximum  depth  of  penetration  in  cm,  m  is  the  mass  in  g,  Ap  is  the 
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presented  area  in  cm2,  V  is  the  striking  velocity  of  a  shard  in  m/s  (Weber  et  al.  2012). 
Based  on  criteria  set  forth  in  Weber  et  al.  2012,  a  fragment  that  penetrates  1  to  12  mm 
causes  mild  injury,  and  a  fragment  that  penetrates  more  than  12  mm  causes  serious  injury. 
Going  forward,  we  must  determine  how  best  to  integrate  the  information  from  more  advanced 
models  into  HENRE  for  casualty  estimation. 

We  have  not  addressed  a  revision  of  blunt  injury  models  here.  However,  work  has  been 
done  correlating  blunt  trauma  to  specific  body  parts  with  injury.  For  example,  White  et 
al.  measured  blunt  trauma  to  the  head  due  to  a  10-lb  missile  as  a  function  of  impact  ve¬ 
locity  (White  et  al.  1965).  Furthermore,  another  group  at  ARA  has  developed  Massive 
Projective/Whole  Body  Displacement  software  which  simulates  collisions  between  the  hu¬ 
man  body  and  large  rigid  objects  to  predict  the  risk  of  injury  due  to  blunt  impact  (Walilko 
et  al.  2014).  This  software  uses  mechanical  principles  to  model  collisions  to  different  parts 
of  the  body.  Future  work  may  involve  incorporating  this  software  into  HENRE. 

In  conclusion,  the  models  presented  here  provide  an  incremental  improvement  to  blast 
injury  calculations  and  facilitate  injury  prediction  in  urban  environments.  By  combining 
the  translation  and  probit  models,  the  probability  of  injury  and  mortality  can  be  predicted 
for  an  individual  exposed  to  complex  nuclear  blast  waves.  As  discussed  above,  a  number  of 
other  improvements  could  be  undertaken  in  the  future. 
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Appendix  A  R  Code 

A.l  Script  for  Running  Models  and  Generating  Output 

The  following  R  script  was  used  to  run  the  models  given  in  Appendix  A. 2  and  A. 3.  The 
script  contains  code  which  verifies  the  R  implementation  and  performs  validations.  This 
code  outputs  the  Figures  given  in  Section  3.  Some  of  these  verifications  and  validations  rely 
on  input  datasets,  including  the  digitized  version  of  Figure  2.2  from  Bowen  et  al.  1961,  Table 
4.1  from  Bowen  et  al.  1961,  Figure  Al  from  Fletcher  et  al.  1975,  Figure  4  from  Otte  2001, 
and  Figure  11.8  from  Needham  2010. 

#Script  for  implementing  translation  models  discussed  in  DTRA -TR - 1 4 -23 
# 

#Jacqui  Wentz 
# jwentz@ara . com 
#1 / 1 4 / 20 1 5 

#  Set  up  - 

setwd("C:/Users/j wentz/Documents/Pro jects/NSF/Injury  Criteria  Modeling/1') 
library (deSolve) 

1 ibra  ry ( ggplot2 ) 
library (scales) 

source("Scripts/objectVelocityBowen  .  R") 

source(" Scripts/ob jectVelocityWithDecelFletcherFinal.R") 

#Conversions 

psi_to_kpa  <-  6.89475729  #Conversion  factor  to  change  from  psi  to  kpa 

#Determine  ratio  of  overpressure  duration  to  dynamic  pressure  duration  (based 
on  Figure  2.2  and  Table  3.1  from  Bowen  1961) 
tratio_data_1  <-  read . csv("Digitized//Bowen1 961 //Bowenl 961 _Figure2point2 . csv”) 
tratio_data_1 $tp_0ver_tu  <-  1 / tratio_data_1 $tu_0ver_tp 
tratio_data_2  <-  data. frame(Ps=c(0.068,  .10, .15, .20,  .25,  .30,  .35  ,  .4, 

.5, .6, .7, .8,1  ,1.3, 1.7), 

tp_0ver_tu  =  c(.9,  .885,  .875,  .855,  .840,  .835,  .805, 

.793,  .76,  .74,  .72,  .71,  .675, .635,  .585)) 
tratio_data_2$tu_0ver_tp  <-  1 / tratio_data_2$tp_0ver_tu 
tratio_data  <-  r bi nd ( t ra t i o_da ta _1 , t ra t io_data _2) 

#Load  in  functions 

inverse  <-  function  (f,  lower=0,  upper=10)  { 

function  (y)  un i root ( ( f unc t i on  (x)  f(x)  -  y) ,  lower  =  lower  ,  upper  =  upper )  [1 ] 

} 

f ind_overpressure  <-  inverse (function  (x)  2 . 5 *x * 2 / ( 7  +  x) * ( 1 +2* 1 0 * ( -8) *x A 4) / 

( 1  +1  0 A ( -8) *x  A  4)  ,  0,  500) 

f ind_decay_rate  <-  i n ve r se ( f unc t ion  (x)  ( exp ( -x ) +x - 1 ) /x A 2 , -1 0 , 1 0) 

#  Verifying  acceleration  portion  of  translation  model  - 

#Parameters  obtained  from  Bowen  1961,  Section  5.2 
ps  <-  5 . 3  #psi 

p 0  <  —  13.3  #psi 

tp_plus  <-  .964  #s 
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c 0  <-  1120  # f t / s 

alpha  <-  0.03  #ft~2/lb 

Ps  <-  ps/p0 

tp_Over_tu  <-  approx (x=tratio_data$Ps , y=tratio_data$tp_Over_tu ,  xout=c (Ps) ) $y 
tu_plus  <-  tp_plus/tp_Over_tu 

A  <-  ( al pha *  1 44) * ( p0* 3 2 . 1  7 4049 ) * tu_pl us / c0  #unitless 
Vini  <-  c(V=0,D=0) 
params  <-  c (Ps , A) 

out  <-  ode(y  =  Vini,  func  =  ob j ect Ve loc i tyBowen  ,  times  =  seq (0 ,  1  , . 002)  ,  parms 
=  params,  method  =  "ode23") 
max (out [  ,  2]) *c0 

#Gives  max  veloctiy  of  23.0  ft/s  compared  with  value  from  paper  of  23.4  ft/s 
#Table  4.1  Verification 

BowenTabData  <-  read.csv("Digitized/Bowen1961/Bowen1961Table4.1  . cs v  ”  ) 

Ps  <-  0.1 

modelData  <-  data  .  frame (P  =  c (), V  =  c (), D  =  c (), A  =  c  ()  ) 
maxVModel  <-  data.frameQ 
Vini  <-  c ( V=0 , D  =  0) 

for  (Ai  in  un i que ( BowenTabDa t a $A) )  { 
params  <-  c ( Ps , Ai ) 

out  <-  ode(y  =  Vini,  func  =  obj ectVeloc i tyBowen  ,  times  =  seq  (0 , 1  ,  .  002)  , 
pa rms  =  pa  rams ,  method  =  "ode23") 
minlndex  <-  min (which (out [,  2]  %in%  max (out [  ,  2] ) ) ) 
out  <-  out [1 : minlndex  ,  ] 

maxVModel  <-  rbind(maxVModel  ,data.frame(P  =  Ps,Vdata  =  max(subset(BowenTabData  ,  A 
==Ai)$V) ,V=out[minIndex ,2] , D=out[minIndex ,3] , A= Ai ) ) 
modelData  <-  rbind(modelData,data.frame(P=Ps,V=out[,2],D=out[,3],A=Ai)) 


ggplot  ( da  ta  =  BowenTabDa ta  )  + 

geom_point(aes(D,V,color=factor(A)),size=2,shape=1)  + 
geom_text (data=maxVModel , 
aes (D , V , 

label  =  pas  te  (  "  A= "  ,  A  ,  "  \  nMax  .  Velocity  (Model):  11  ,  round  (V  ,  5)  , 

”\nMax.  Velocity  (Data):  ",round(Vdata,5),sep=”"), 
color=factor(A))  , 
v j us t  =  .  5 , h j us t  =  -0 . 06 , size  =  3)  + 
geom_line(data=modelData , aes(D,V, color=factor(A)))  + 
scale_x_continuous(limits=c(0,0.02))  + 
scale_y_continuous(limits=c(0,0.05))  + 
xlab(” Distance  (dimensionless)”)  + 
ylab(” Velocity  (dimensionless)”)  + 
scale_colou r _br ewe r ( pa le t te  =  ” Set  1  ")  + 
theme_bw()  + 

theme(text=element_text(size=12), 
legend . position  =  "none”  , 

panel . grid . major=element_line(color="gray") , 
panel . grid . minor=element_blank()) 
#ggsave("Figures/BowenVerification_v7.pdf”,width=6.5,height=5.5) 
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#  Verifying  Fletcher  Report  Run  10,  10kT 


#Baseline  parameters 
p0  <-  1 4 . 7 * ps i _ to_kpa  #kpa 
c0  <-  340 .29  #m/s 
T0  <-  298  #K 

R  <-  287.058  # J/ kg/K  =  mA2/s“2/K 
g  <-  7/5  #1;  adiabatic  index 
m0  <-  75  #kg 

#Blast  wave/object  parameters 


W  <-  10  #kT  ; 

qs  <-  1 4*psi_to_kpa  #kPa 

ps  <-  27 . 4*psi_to_kpa  #kPa 

tp_plus  <-  0 . 1 65*W A  (1 /3)  #s; 

tu_plus  <-  0 . 3 1 2*W * (1 / 3)  #s; 


Ip  <-  1 . 1 5 *W A ( 1 / 3) * ps i _to_kpa 
Iq  <-  0 . 474*W A ( 1 / 3) * ps i _ to_kpa 
m  <-  75  #kg 


Given  in  Fletcher  report 
;  Given  in  Fletcher  report 
;  Given  in  Fletcher  report 
Given  in  Fletcher  report 
Given  in  Fletcher  report 
#kPa  s;  Given  in  Fletcher  report 
#kPa  s;  Given  in  Fletcher  report 


#Dimensionless  quantities 
RBar  <-  R*T0/c0'2 
Ps  <-  ps/p0 
Qs  <-  qs/p0 
M  <-  m/m0 

#Acceleration  coefficient  for  man  initially  prone  with  a  random  orientation 
Dlim  <-  1 . 79832/tu_plus/c0  #Distance  at  which  acceleration  coefficient  is  Amax 
Amin  <-  . 0028 6 74260 1  * ( p0* 1  000) * tu_pl us / c0 
Amax  <-  . 006 1 444843 1 *( p0* 1 000) * tu_pl us / c0 

Achange  <-  tu_pl us A 3 *c0* ( p0 *  1  000) / 3 1  . 65 A 2  #Determines  how  A  changes  with  DA2 


#Deceleration 

FI =8 . 9029* ( . 3048/ c0) A . 61 692*tu_plus 


#Find  wave  decay  rates  based  on  impulse  and  peak 
n  <-  f i nd_decay _ ra te ( I p / p0 / tp_plus / Ps ) $  root 
r  <-  f i nd_decay_ rate ( Iq / p0 / tu_plus / Qs )$ root 


#Run  ODE 

Vini  <-  c ( V=0 , D=0) 

params  <-  c(Ps,Qs,n,r,tp_plus/tu_plus,M, RBar, Dlim, Am in, Amax, Achange, FI ,g) 
out  <-  ode(y  =  Vini,  func  =  objectVelocityWithDecelFletcherFinal  , 
times  =  seq(0,5,.002),  parms  =  params,  method="ode23”) 
modelData  <-  data . frame(time=out[ , 1]*tu_plus , velocity=out[  , 2] *c0  , 

distance=out[ ,3]*c0*tu_plus) 


#Load  in  and  set  up  verification  data 

Curve  <-  read  .  csv(''Digitized/Fletcher1  975/FigureAI  Curve  .  csv") 

TimePoints  < -  read . csv("Digitized//Fletcher1 975//FigureA1TimePoints . csv") 
TimePoints$label2  <-  paste(factor(sprintf("%.2f",round(TimePoints$label  ,2)))," 
s”  ,  s  e  p  =  "  ") 

modelTimePoints  <-  data. frame(time=c() ,velcoity=c() ,distance=c()) 
for  (time  in  T i mePo i n ts $  label )  { 

modelTimePoints  <-  rbind (modelTimePoints , modelData [which (abs (modelData$time - 
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time)==min(abs(modelData$time-time))),]) 

} 

#Generate  comparison  plot  between  data  and  simulation 
plot  <-  ggplotQ  + 

geom_line(data=Curve , aes(Distance.m, Velocity. mpers))  + 
geom_point(data=TimePoints,aes(Distance.m, Velocity. mpers), shape=1)  + 
geom_text(data=subset(TimePoints , label2!="0.20  s") , 

aes (Distance. m, Velocity. mpers , label=label2)  , 
v j us t  =  . 1  ,  h j us t  =  - . 2 , s i ze  =3 )  + 
geom_text(data  =  subset(TimePoints , label2==”0.20  s")  , 

aes (Distance. m, Velocity. mpers , label  =  label2)  , 
v j us t  =  . 2 , h j us t  =  1 . 2 , s  i  ze  =3  )  + 

geom_line(data=modelData,aes(distance, velocity), color=" red", s i ze  =  . 25)  + 
geom_point(data=modelTimePoints, 

aes (distance , velocity) , 
color="red",shape=1)  + 
scale_x_continuous(limits=c(0 ,23))  + 
xlab(" Distance  Traveled  (m)")  + 
ylab(" Velocity  (m/s)”)  + 
theme_bw()  + 

theme(text=element_text(size=12), 

panel  .  grid . major  =  element_line(color  =  ”gray") , 
panel . grid . minor=element_blank()) 
mergedTimePoints  <-  cbind (modelTimePoints , TimePoints) 
for  (i  in  1 : n row ( me rgedT i mePo i n t s ) )  { 
row  <-  mergedTimePoints  [i  ,  ] 
plot  <-  plot  + 

geom_line(data=data.frame(velocity=c(row$velocity,row$Velocity. mpers), 

distance  =  c(row$distance , row$Distance.m))  , 
aes (distance , velocity) , linetype=2) 

} 

plot 

#ggsave("Figures/FletcherVerification_v7.pdf",width=6.5,height=4.5) 

#  Model  validation  using  accident  data  from  Otte2001  - 

#Simulation  parameters 
c0  <-  340.29  #m/s 

tu_plus  <-  1  #Need  to  have  value  as  place  holder,  doesn’t  affect  results 
FI  <-  8 . 9029* ( . 3048/C0) A . 61 692*tu_plus 

#Run  simulation  for  different  initial  velocities  and  store  max  distance 
decelData  <-  da ta . frame ( model  =  c (), ca lc  =  c  () ) 
for  (vi  in  seq(1  ,30,1))  { 

V  <-  vi/c0 

Vini  <-  c ( V=V , D=0) 

params  <-  c (0 , 0 , 1  , 1  , 1  , 1  , RBar  ,0 , 0 , 0 , 0 , FI  ,  g)  #Set  blast  wave  params  to  zero 
out  <-  ode(y  =  Vini,  func  =  ob j ec t Ve loc i ty Wi t hDece 1 F le t che r F i na 1  , 
times  =  seq(0,10,.002),  parms  =  params,  method=”ode23") 
modelData  <-  data . frame(time  =  out[ , 1]*tu_plus , velocity  =  out[  ,2]*c0 , 

distance=out[ ,3]*c0*tu_plus) 

decelData  <-  rbind(decelData ,data. frame(vi=vi , model=max(modelData$distance) , 

calc  =  1 0‘ ( -1  .15822  +  1  . 61692*1 ogl 0 (vi)) 
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} 


)) 


#Load  in  validation  data 

otteData  <-  read.csv("Digitized//Otte2001//Figure4.csv") 

#Plot  data  against  model  predictions 
ggplot(data=decelData ,aes(vi , model))  + 
geom_line  ()  + 

geom_point(data  =  otteData , aes(x*.28 , Curvel ) , color  =  ”blue"  ,  shape  =  1 )  + 
xlab('' Velocity  (m/s)")  + 
ylab(" Distance  (m)")  + 
theme_bw()  + 

theme(text=element_text(size=12), 
legend . position  =  ”none"  , 

panel . grid . major=element_line(color="gray") , 
panel . grid . minor=element_blank()) 

#ggsave("Figures/0tteDecelerationValidation_v7.pdf”,width=6.5,height=4.5) 

#  Model  Validation:  Jeep  Data  ( Needham20 1 0 )  - 

#Load  in  validation  data 

jeepData  < -  read . csv(”Digitized/Needham201 0//Figure1 1 point8 . csv”) 

#Baseline  parameters 
T0  <-  298  #K 

R  <-  287.058  # J/ kg/K  =  mA2/s“2/K 
RBar  <-  R*T0/c0“2 
g  <-  7/5 

p0  <-  1 4 . 7*psi_to_kpa  #kPa 
c0  <-  340.29  #m/s 
m0  <-  75 

#Object  parameters 

alpha  <-  .005  #m"2/kg;  very  rough  approximation  for  a  jeep 
m  <-  1000  #kg;  approximate  weight  of  jeep 

#Dimensionless  parameters  (same  for  each  simulation) 

M  <-  m/m0 

Dlim  <-  -1  #Set  Dlim  to  less  than  one  so  that  A  is  constant  at  Amax 
Amin  <-  0 
Achange  <-  0 


#Iterate  through  dynamic  impulses  and  run  simulations  under  several  conditions 
.  Store  max  and  min  for  each  dynamic  impulse. 
jeepModelData  <-  data . frame (Iq  =  c () , dmin=c () , dmax=c () ) 

Vini  <-  c(V=0,D=0)  #initial  conditions  of  ode 

for  (Iq  in  c  (  .  5  ,  .  6  ,  .  7  ,  .  8  ,  1  ,  1  .  5  ,  2  ,  2  .  75  ,  4  ,  5  ,  8  ,  1  0 , 1  5  ,  20  ,  27  .  5  ,  50 , 80 , 1  00)  )  { 
dmax  <-  0 
dmin  <-  99999 
print (Iq) 

for  (tu_plus  in  c ( . 1  ,  . 3  ,  .  5  ,  .  7  ,  .  9  ,  1  )  )  { 

FI  <-  8.9029*(.3048/c0) n .61692*tu_plus 
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A  <-  alpha*p0*1 000*tu_plus/c0 
for  (r  in  c  (  .  1  ,  .  4 ,  .  7  , 1  , 3 , 5)  )  { 

Qs  <-  ( Iq / p0 / tu_plus ) / ( exp ( - r ) +r -1 ) * r  "  2 
Ps  <-  f ind_overpressure (Qs) $root 
for  (tp_plus  in  c ( . 1  ,  . 3  ,  .  5  ,  .  7  ,  .  9  , 1 )  )  { 
for  (n  in  c ( . 1  ,  .  4  ,  .  7  , 1  , 3 , 5) )  { 

params  <-  c(Ps,Qs,n,r,tp_plus/tu_plus,M,RBar,Dlim,Amin,A,Achange,F1  , 
g) 

out  <-  ode(y  =  Vini  ,  func  =  ob j ec t Ve loc i ty Wi t hDece 1 F 1 e t che rF i na 1  , 
times  =  seq(0,10,.002),  parms  =  params,  method="ode23") 
dmin  <-  min (dmin , max (out [ , 3] *tu_plus*c0) ) 
dmax  <-  max (dmax , max (out [ , 3] *tu_plus*c0) ) 

} 

} 

} 

} 

jeepModelData  <-  rbind(jeepModelData,data.frame(Iq=Iq,dmin=dmin,dmax=dmax)) 


#Plot  simulation  prediction  against  validation  data 
ggplot ( da ta  =  j eepMode IDa t a  )  + 

geom_ribbon (aes(Iq , dmin , ymin=dmin , ymax=dmax) , a lpha  = . 5)  + 

geom_point(data=jeepData , aes(x*6. 89476, Curvel*. 3048) ,color="red”)  + 

scale_x_log10(breaks=c(1  ,10,100),limits  =  c(.5,100))  + 

scale_y_log1 0 (breaks  =  c  (1  ,10,1  00,1000))  + 

xlab (" Dynamic  Pressure  Impulse  (kPa  s)")  + 

ylab('' Displacement  (m)")  + 

theme_bw()  + 

theme (panel . grid . major  =  element_line (color  =  "gray")  , 
panel . grid . minor=element_blank()) 
ggsave("Figures/jeepDataValidation_v7.pdf",width=4.5,height=4.5) 

#write . csv(modelData , "jeepValidationData . csv") 

# jeepModelData  <-  read . csv (" jeepValidationData . csv  ") 
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A. 2  Bowen  1961  Translation  Model 


The  following  R  function  is  used  to  calculate  an  object’s  maximum  velocity  based  on  the 
model  given  in  Bowen  et  al.  1961.  This  function  is  called  using  an  ODE  solver. 

ob jectVelocityBowen  <-  function(Z,  y,  p)  { 

’Function  calculates  the  distance  and  velocity  derivatives  for  running  the 
translation  model  (called  with  an  ODE  solver).  Based  on  the  Bowen  1961 
translation  model . 

Inputs  : 

Ps  -  peak  overpressure  (dimensionless) 

A  -  acceleration  coefficient  (dimensionless) 

#State  variables 
V  <-  y  [1 ] 

D  <-  y [2] 

#Parameters 
Ps  <-  p  [1 ] 

A  <-  p  [2] 

if  (Ps  <  0.6)  { 

J  <-  1.1 86  *  Ps A (1 /3) 

}  else  if  (Ps  <=  1 )  { 

J  <-  1 
}  else  { 

J  <-  10A4*PsA(-1 / 4)  /  ( 1  0  A  4  +  Ps  A  2) 

} 

K  <-  1-J 

gamma  <-  1/4  +  3.6*PsA(1/2) 

delta  <-  7  +  8*PsA(1/2)  +  2*PsA2/(240  +  Ps) 

Qs  <-  2.5*PsA2/(7+Ps)*(1+2*10A(-8)*PsA4)/(1+10A(-8)*PsA4) 

Us  <-  Ps/ (1 +Ps A (1 /2) ) 

Xdot _Max  <-  3/5*Us  +  ( 1 + ( 3 / 5*Us ) A 2) A ( 1 / 2) 
nu  <-  PsA(1 /3)+.0032*PsA(3/2) 

U  <-  Us*  ( 1 ~(Xdot_Max*Z-D)/Xdot_Max)*exp(-nu*(Xdot_Max*Z-D)/Xdot_Max) 

#U  <-  Us* ( 1  - Z) *exp ( - nu *Z) 
if  (V>=U)  { 

return (list(c(V=0 , D=0 , Zprime=0) )) 

}  else  { 

Q  <-  Qs*(1 -(Xdot_Max*Z-D)/Xdot_Max)*(J*exp(-gamma*(Xdot_Max*Z-D)/Xdot_Max) 
+K*exp(-delta*(Xdot_Max*Z-D)/Xdot_Max)) 

#Q  <-  Qs * (1  - Z) * (J*exp ( -gamma *Z) +K*exp (- del ta *Z) ) 
dV  <-  Q*A*((U-V)/U) A2 
dD  <-  V 

return(list(c(V=dV,D=dD))) 

} 

} 
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A. 3  Fletcher  1966/1975  Translation  Model 

The  following  R  function  is  used  to  calculate  an  object’s  velocity  and  distance  traveled  based 
on  both  acceleration  due  to  the  blast  wave  and  deceleration.  This  function  is  called  using 
an  ODE  solver. 

ob j ec t Ve loc i ty Wi t hDece 1 F 1 e t che r F i na 1  <-  function (Z ,  y,  p)  { 

’Function  calculates  the  distance  and  velocity  derivatives  for  running  the 
translation  model  (called  with  an  ODE  solver). 

Based  on  Fletcher  1966/1975’ 

#State  variables 
V  <-  y  [1 ] 

D  <-  y  [2] 

#Parameters 
Ps  <-  p  [1 ] 

Qs  <-  p  [2] 
n  <-  p  [3] 
r  <-  p [4] 
tratio  <-  p [ 5 ] 

M  <-  p  [ 6 ] 

RBar  <-  p  [ 7 ] 

Dlim  <-  p  [  8 ] 

Amin  <-  p  [ 9 ] 

Amax  <-  p  [  1 0 ] 

Achange  <-  p  [  1 1  ] 

FI  <-  p [ 1 2] 
g  <-  P [1 3] 

#Acceleration  coefficient 
if  (D<Dlim){ 

A  <-  Ami n + Achange *D A 2 
}  else  { 

A  <-  Amax 

} 

#Pressure  calcuation 
Q=Qs*(1 -Z)*exp(-r*Z) 

P=Ps*(1-Z/tratio)*exp(-n*Z/tratio) 
if  (Q<0)  Q  =  0 
if  (P<0)  P=0 

#Calculating  wind  velocity  within  shock  wave  using  adiabatic  conditions 
Ts  <-  ( 7+ Ps ) / ( 7+6 *Ps ) * ( 1 + Ps )  #Using  RH  relation 
TBar  <-  Ts* ( ( P +1 ) / ( Ps +1 ) ) A ( ( g - 1 ) /g) 

RhoBar  <-  (P+1)/TBar  #RhoBar  =  rho  *  RT0/P0 
U  <-  sqrt (2*Q*RBar/RhoBar ) 
if  (is.na(U))  U  =  0 

#Return  changes  to  velocity  and  distance  traveled 
if  (V<0)  { 

return(list(c(0,0))) 

}  else  if  (V>=U)  { 
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dV  <-  -1 *F1 *(V*(1 /M) A (1 /6)) A0. 38308 
dD  <-  V 
}  else  { 

dV  <-  Q*A* ((U-V)/U)  "2  -  FI *(V*(1 /M) A(1 /6)) ‘0. 38308 
dD  <-  V 

} 

return (list(c(dV , dD))) 
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Abbreviations,  Acronyms,  and  Symbols 


AIS 

AmedP-8 

ARA 

Cl 

DT 

DTRA 
EM-1  Ch.  14 
FI 

FME 

HENRE 

HPAC 

HRP 

IPI 

ISS 

LFMER 

SS 

MHGP 

MI 

NATO 

ORCA 

PDCALC 

PI 

PRCC 

R 

SFOM 

SI 


abbreviated  injury  scale 
Allied  Medical  Publication  8 
Applied  Research  Associates,  Inc. 
combat  ineffectiveness 
decelerative  tumbling 
Defense  Threat  Reduction  Agency 
Effects  Manual- 1  Chapter  14 
fatal  injury 

flexible  modeling  environment  for  modeling,  sensitivity,  and  Monte  Carlo 
analysis  (R  plug-in) 

Health  Effects  from  Nuclear  and  Radiation  Environments 

Hazard  Prediction  and  Assessment  Capability 

Human  Response  Panel 

immediate  permanent  ineffectiveness 

Injury  Severity  Score 

Lovelace  Foundation  for  Medical  Education  and  Research 

Severities  Sum 

Multi-Hit  Glass  Penetration 

moderate  injury 

North  Atlantic  Treaty  Organization 

Operational  Requirement-based  Casualty  Assessment 

Probability  of  Damage  Calculator 

perpendicular  impact 

Personnel  Risk  and  Casualty  Criteria 

software  programming  language  for  statistical  computing 

Shard  Fly-Out  Model 


serious  injury 


58 


DISTRIBUTION  LIST 
DTRA-TR-1 5-23 


DEPARTMENT  OF  DEFENSE 


DEFENSE  THREAT  REDUCTION 
AGENCY 

8725  JOHN  J.  KINGMAN  ROAD 
STOP  6201 

FORT  BELVOIR,  VA  22060 
ATTN:  P.  BLAKE 

DEFENSE  TECHNICAL 
INFORMATION  CENTER 
8725  JOHN  J.  KINGMAN  ROAD, 
SUITE  0944 

FT.  BELVOIR,  VA  22060-6201 
ATTN:  DTIC/OCA 

DEPARTMENT  OF  DEFENSE 
CONTRACTORS 

QUANTERION  SOLUTIONS,  INC. 
1680  TEXAS  STREET,  SE 
KIRTLAND  AFB,  NM  87117-5669 
ATTN:  DTRIAC 


DL-1 


